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ith (current or updated) average 
ith sample value 

temperature-shift correction constant 

corrected bridge output voltage 

measured bridge output voltage 

fictive no-flow bridge voltage 

average value of the measured bridge voltage 

sample number 

constant 

number of samples 

log slope (about 0.45, typical value) 

overheat ratio (usually about 1.8) 

average of quantity a, defined by equations (15) 

sum of a(i) from i = 1 to i = Ns 

current flow temperature 

calibration flow temperature 

calibration velocity 

effective cooling velocity 

the vector of effective cooling velocities 

instantaneous velocity components in x,y,z directions 

instantaneous velocity vector 

magnitude of the velocity vector 

Cartesian coordinate directions 

wire resistivity coefficient (typical value for tungsten = 0.005 per °C) 
yaw calibration angle setting 
yaw angle 
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[ip ] angular sensitivity matrix 

Subscripts: 

1,2 used to denote a particular quantity referred to wire 1 or wire 2 

u v , uw quantity determined in uv or uw measurement plane 

Superscripts: 

1 fluctuating quantity, e,g. , u = u + u f , or used to indicate an inter- 

mediate quantity (e.g., eq . (11)) 

time average 
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SUMMARY 


The report describes a system for rapid computerized calibration acquisition, 
and processing of data from a crossed hot-wire anemometer. Advantages of the system 
are its speed, minimal use of analog electronics, and improved accuracy of the result- 
ing data. Two components of mean velocity and turbulence statistics up to third 
order are provided by the data reduction. The report presents details of the hardware, 
calibration procedures, response equations, software, and sample results from measure- 
ments in a turbulent plane mixing layer. 


1 . INTRODUCTION 


A system for rapid computerized acquisition and processing of analog signals 
from a hot-wire anemometer is described. Most of the discussion considers the 
crossed-wire anemometer; the simpler case of a single wire is briefly addressed. The 
objective of the work was to develop a system for measuring the statistical properties 
of the velocity field of a moderately turbulent, two-dimensional, incompressible, 
isothermal air flow. This report is intended to document the system for reference in 
several researches in which the system will be applied. 

A crossed hot-wire probe operated by constant-temperature bridges (CTA) can 
provide a signal from which the two instantaneous components of velocity and related 
statistical properties of a moderately turbulent flow may be derived. Traditionally, 
analog hardware is used to linearize and process the signals from multichannel CTA 
systems. More analog hardware is required to proces s and av era ge th e resulting sig- 
nals to obtain the desired statistics; for example, u ?2 ,v 12 and u f v f for a crossed- 
wire. The present report describes a more modern approach wherein all nonlinear ana- 
log processing hardware is replaced by a computerized system for probe calibration, 
data acquisition, and reduction. The approach described differs from the old fully 
analog methods in that it uses fewer, more stable electronic components, and, more 
importantly, is much faster. Complete calibration and acquisition for a typical 
20-point profile required less than 60 min. Corrections (e.g., of ambient temperature 
drift) are easily implemented in software. Although specific hardware with particular 
response equations are presented, the approach taken is very general and could be 
easily adapted to different measurement conditions, other types of anemometers, or 
three-wire probes. 

The experimental procedure and apparatus are discussed in detail in the next 
section. Hot-wire response equations — which give simple but accurate relations for 
the cooling law, directional sensitivity, and the effects of ambient temperature 
drift — and algorithms for computing desired statistical signal properties, such as 
mean, variance, and higher-order cross-correlations, are presented in the third 
section. Sample results from an experiment to measure turbulence quantities in a 
plane mixing layer are provided in the fourth section, and concluding remarks are 
presented in the final section. Software for both single- and crossed-wire systems 
written in BASIC to run on the HP 9845B desk-top computer is included in the appendix. 



2. EXPERIMENTAL PROCEDURE AND APPARATUS 


Experimental Procedure 

The hardware configuration of the system is conceptually simple and requires 
little routine adjustment. Two DISA constant-temperature bridges were employed 
to drive the crossed-wires . A fixed dc shift and gain were applied to the bridge 
outputs which were low-pass filtered (to remove high-frequency electrical noise) 
and then input to a bipolar 12-bit A/D converter. An external clock generated 
a sample pulse (typical sampling rate was 500 Hz) to the A/D which caused the two 
inputs to be frozen ("sample-and-hold") . These channels were sequentially converted; 
the converted values passed through a multiplexer to a high-speed, 16-bit parallel 
interface to the computer. After filling the computer’s data buffer, raw data were 
written to floppy disk. The raw data were reduced off-line to provide the signal 
statistics. Before data acquisition began, a complete system calibration check was 
performed, and relevant quantities were written to the floppy-disk data files. Taking 
up to 3,000 samples per data point required less than 2 min; the calibration procedure 
required about 15 min. Thus, a typical 20-point profile required less than a hour for 
data acquisition (including calibration) , using the program shown here (see appendix 
sec. A. 1.2). The current off-line data reduction program (see appendix sec. A. 1.3) 
took about 45 min to reduce these data. 

Simple but established models have been implemented to describe the sensitivity 
of the constant-temperature hot-wire to variations in flow velocity, wire orienta- 
tion, and ambient temperature drift. The relations selected are applicable to 
incompressible, isothermal flow of "clean" (filtered) air over a fine wire with 
"moderate" local turbulence intensity (less than about 30%) and negligible instan- 
taneous local flow reversal frequency. King’s law was used to provide a relation 
between the "effective" cooling velocity and the bridge output. A "cosine law" was 
used to relate the effective cooling velocity to the magnitude and direction of the 
velocity vector. The entire calibration was shifted to account for the effects of 
small ambient temperature changes over the course of a run, based on the overheat 
ratio and wire resistivity. Once the King’s law calibration was implemented and 
shifted for ambient temperature drift, the effective velocities measured by each wire 
were used to solve for the instantaneous velocity components in the measurement plane 
by inverting the angular sensitivity matrix. Then the instantaneous values of the 
velocity components were used to form the statistical properties of the signal. Com- 
puted results included the_mean and vari a nce, as wel l as se cond- and third-order 
cross-correlations (i.e., u, v, u’ 2 , v’ 2 , u’v’ , u’ 2 v, v’ 2 u’ for wires in the u-v 
plane) . 


Experimental Apparatus 

Figure 1 shows the system hardware schematically and provides a list of compo- 
nent manufacturers and model numbers. In the following discussion of the hardware, 
the system is divided into four areas: (1) probes used, (2) DISA bridges and 

signal conditioners, (3) NASA LDV-A/D computer interface, and (4) HP computer inter- 
face, desk-top computer, and floppy-disk drive. Little explicit mention will be made 
of the digital voltmeter or oscilloscope , which are used to monitor the analog sig- 
nals, or of the pulse generator, which provides a sample trigger to the A/D converter. 

Crossed-wire probes- Miniature crossed-wire probes were manufactured in-house at 
Ames. These probes had two nominally perpendicular wires mounted at angles of about 


2 



±45° to the probe stem (see fig. 2). The planes that contain each wire and its two 
support needles were parallel and separated by about 1 mm. The wires were 5-ym 
tungsten elements about 1 mm long; they were welded to the supports. The "cold" 
(room-temperature) resistance of each wire was about 5 ohms, and the operating over- 
heat ratio used was 1.8 (ratio of hot, or operating, resistance to cold, or unheated, 
wire resistance). The probe tips were mounted in a stem which was in turn held in a 
rotating collar. A spring-loaded mechanism in the collar could be rotated to yaw the 
probe in five fixed positions in 5 0 increments. In this way, yaw calibration of the 
wire was performed. Often it was convenient to mount the probe so that the axis of 
rotation of the collar was parallel (rather than perpendicular) to the desired plane 
of measurement. Then a special procedure was used wherein the yaw calibration was 
performed in the u-w plane, but the probe was then rotated about the stem axis for 
measurement in the u-v plane (see Summary: Calibration and Dat a-Acquisition Proce- 
dure in the next section). 

Bridges and signal conditioners - Two DISA 55D31 bridges were used to operate the 
crossed hot-wires at constant temperature. The cable resistance was compensated for 
using a shorting probe. The probe resistance was measured, and the operating resis- 
tance of each probe was set at 1.8 times the cold resistance; this is a nominal value, 
because the operating resistance was not changed to account for daily changes in 
ambient temperature. However, the frequency response was checked daily as described 
below. 

The frequency response of each channel was optimized using the internal 1-kHz 
square-wave generator and varying the bridge gain (as well as the cable compensation 
adjustments). Typically, the "3-dB-down" point was about 20 kHz on each channel 
(Freymuth method). A bridge gain of 4 and an HF filter setting of 2 were usual values. 
The bridge output voltage ranged from about 3 V (no-flow) to about 4.7 V (25 m/sec 
flow velocity); since the input voltage range of the available A/D converter was 
-10 to +10 V (with 12-bit resolution), linear signal conditioning elements were used 
to obtain better resolution of the bridge output voltage. 

Commercial signal conditions manufactured by DISA (model 55D26) were used to dc 
shift, to amplify, and to low-pass filter the bridge outputs before A/D conversion. 

A dc shift of about 4 V and a gain of 10 were applied to each channel, yielding a 
voltage covering most of the ±10-V range of the 12-bit A/D converter. Note that the 
exact values of the gain and offset used were calibrated during setup of the system 
for each run (see the description of the software below), so that it was not necessary 
for the actual gain or offset to be determined from the nominal settings. The low- 
pass filter cutoff was 10 kHz with 18 dB/octave roll-off; the high-pass stage was set 
to "direct." The filtering was intended to eliminate spurious electrical noise which 
would contaminate measurements of low turbulence intensity. 

NASA LDV-A/D computer interface- Two channels of the A/D conversion stage of the 
NASA LDV-A/D computer interface (ref. 1) were used to convert the conditioned signals 
to digital form and multiplex them to the HP computer interface (described below). 
Figure 1 shows how the NASA LDV-A/D fits into the measurement system, and figure 3 
gives the details of the various settings and connections for the interface. A sample 
pulse was generated by a Tektronix PG 508 pulse generator and applied to the "chan- 
nel 1 event" and "external reset" inputs. This signal acted as a sampling trigger to 
initiate a "sample- and-hold" of the two A/D inputs when the interface was enabled by 
the computer setting "CTL0." 

Four 16-bit TTL data words were then multiplexed to the computer. The first word 
gave the time between samples (a count dependent on the operating mode (see ref. 1 for 
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details)); these data were typically discarded after verifying that the nominal sample 
rates were actually attained. The second word gave status information; it was also 
checked then discarded. The third and fourth words were the two channels of anemom- 
etry data in directly computer-compatible integer representation. Since the A/D that 
was used had only 12 bits of resolution, this required that the upper 4 bits be either 
all 0’s or all l’s, depending on whether the converted voltage was positive or nega- 
tive, respectively. Thus, the converted A/D data could take on integer values of 
-2,048 to +2,047, inclusive. 

An important feature of the NASA LDV-A/D hardware was the fast sample-and-hold 
hardware. The measurement strategy required that the two channels be sampled simul- 
taneously so that the instantaneous-velocity vector components in the measurement 
plane could be determined. Since the frequency response of the anemometer was limited 
to about 20 kHz, it was desirable that the two signals be sampled within, say, 10 jisec 
or so. The sample-and-hold hardware in the LDV-A/D was capable of locking in the 
input analog signal within 0.5 psec of receiving the sample trigger. Then the A/D 
conversion and multiplexing could take place asynchronously (as long as all four words 
were passed before the next sample pulse). For each A/D conversion, a minimum of 
14 psec were required. 

HP computer and parallel interface- Multiplexed data were passed to the HP com- 
puter from the NASA LDV-A/D using the HP 98032A high-speed 16-bit parallel interface. 
Jumpers labeled "9,B,D M were connected inside the 98032A for proper operation with 
the LDV-A/D. A select code of 10 (screw setting on the 98032A) was set for use with 
the software described in the appendix. 

A data buffer of 24 kbytes was provided in the memory of the HP 9845B desk-top 
computer for storage of up to 3,000 samples obtained from one measurement location 
(four words of 2 bytes each are passed from the LDV-A/D). Since two of the four data 
words passed for each sample were merely monitored and discarded as described above, 

12 kbytes of raw data remained to be stored for each point. An HP 9895A floppy-disk 
drive was used for archival storage of raw data. The buffered data were written in 
real time to a sequential-access floppy-disk file. Enough header information was 
written to each file to identify the run, as well as to reproduce calibration tables 
and correct for ambient temperature drift. About 1 min per point was required for 
storing the data on floppy disk. 

Fairly rapid and simple data buffering was possible with the HP computer, using 
convenient high-level commands. Sampling rates as high as 10,000 samples per second 
were used (indicating interface transfer rates 4 times as high in 16-bit words per 
second). Direct memory access (DMA) was not used; the processor was simply dedicated 
to the real-time task. 

The system described above was also streamlined for the simplified case of single- 
channel, hot-wire anemometry. Of course, just half the hardware shown before the A/D 
in figure 1 was needed. The calibration procedure consisted of simply setting up the 
signal conditioner and compiling the King’s law data. On-line data reduction was 
implemented, since this could be done rapidly with only one channel. Thus, the 
look-up table was constructed immediately following the static velocity calibration 
and implemented after filling each data buffer. Mean and fluctuating values of 
velocity were stored; raw data were not usually archived. 

A D1SA 55P11, platinum-plated, tungsten hot-wire was used for the single-wire 
work. The overheat ratio and bridge setup were exactly as described earlier. The 
BASIC program to run the single-wire anemometry system is provided in the appendix. 


4 



3. CROSSED-HOT-WIRE RESPONSE EQUATIONS 


Simple relations that describe the response of a fine, heated wire to varia- 
tions in flow velocity, orientation, and ambient temperature drift have been incor- 
porated into the off-line data reduction program. Of course, the equations selected 
have a strong effect on the type of calibration performed and on the accuracy of the 
results. The response equations are first discussed separately, the calibration 
procedure is summarized, and the algorithms actually used for implementing the 
response equations and computing the signal statistics are then given. 


Static Velocity Response: King’s Law 


The effective, instantaneous wire velocity was assumed to be related to the 
bridge output voltage by the generalized King’s law: 


U „ = K(E%. - E z ) 
eff eff o 


2 x i/n 


( 1 ) 


This is an approximate relation which has been determined to be fairly accurate in 
describing the steady-flow heat loss over cylinders in cross-flow; the (constant) 
value of the log slope n selected (0.45) provides a good fit to experimental data 
at moderate-to-low Reynolds numbers, based on wire diameter. Our steady-flow calibra- 
tion data fit King’s law with an rms error of 0.5% over the range of 5-25 m/sec. 
Although simple interpolation or polynomial fit of the calibration may seem equally 
acceptable, the log-linear King’s law fit provides a quick way to evaluate whether 
the calibration is ’’typical,” and it smooths minor ’’jitter” in the calibration data. 
Also, the King’s law calibration may be confidently extrapolated slightly outside the 
range of the actual calibration data. 


Yaw Sensitivity: Cosine Law 

The ’’effective cooling velocity” is taken to be that component of the velocity 
vector perpendicular to the wire. This assumption implies that a wire yawed in a 
uniform constant -velocity stream will respond to an effective flow velocity that is 
proportional to the cosine of the yaw angle (see fig. 2 for nomenclature). Neglecting 
the axial component of velocity (along the wire axis) is strictly an approximation — 
one that is often made, however, and one that works quite well for moderately turbu- 
lent flows. 


U eff = | V | cos ip (2) 

For two wires at angles \p 1 and ip 2 to the reference coordinate direction, the follow- 
ing two equations result, if sensitivity to out-of-plane velocity fluctuations is 
neglected: 

U effl “ u cos ip i + V sin ijj 1 = Kj (E* - E^) l/ni 

/ ( 3 ) 

u e ff 2 = U cos <J> 2 + v sin \p 2 = K 2 (E% - Eq 2 ) 1/n 2 

Here u and v are the instantaneous velocity components in the measurement plane. 

Note that U e £f for each wire may be found immediately from the King’s law 


5 



calibration; only the bridge voltage must be known. A matrix form more convenient for 
discussion or generalization is 


[u eff ] = M[V] 


(4) 


The entries for the angular sensitivity matrix are found by a yaw-calibration proce- 
dure described below. Then, the inverse of this matrix is computed once and stored. 
Solution for the instantaneous velocity components is simple: 

[V] = [^]" 1 [U eff ] (5) 


The explicit solution of equation (3) gives the equations actually used to compute the 
instantaneous velocity components: 


(sin ^)U eff2 - (sin ^ 2 )U effi 
(sin 4 ^) (cos \p 2 ) - (cos 4 ^) (sin \p 2 ) 

> 

(cos ^ 2 )U effl - (cos ^ x )U eff2 
(sin ^ 1 )(cos \p 2 ) - (cos i^Msin \p 2 ) 


( 6 ) 


Temperature Drift Correction 

Temperature drift of a few degrees Fahrenheit is commonly encountered and has 
some effect on measurement accuracy since the probe is operated at constant tempera- 
ture. For typical wire-overheat ratios and sensor resistivity properties, the tem- 
perature difference between the wire and flow is a few hundred degrees Fahrenheit, so 
that a few degrees drift will change the perceived heat-transfer coefficient between 
wire and flow by a percent or so. A small correction to the bridge voltage is 
applied to account for this variation in ambient temperature: 


E - E 
corr meas 

E 

meas 


2 (OHR) 


(T - 


T -) E 
cal 7 


(7) 


where C represents the percent shift required for the instantaneous bridge voltage — 
fixed by the operating parameters and current-flow temperature. Since fluctuations in 
the bridge voltage are fairly small, a further simplifying approximation is that the 
instantaneous bridge voltage can be simply shifted by the fixed percentage of the 
average output voltage. Then the calibration is easily implemented after the shift is 
computed and applied to each voltage reading. The equation actually used is then 

E = E + C x E (8) 

corr meas meas 


Note that each voltage reading is corrected for ambient temperature changes. The 
correction will, therefore, influence both mean and fluctuating time-averaged results, 
as it should. 
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Summary: Calibration and Data-Acquisition Procedure 

Figure 4 shows a flowchart that represents the verbal description of the calibra- 
tion and data-acquisition procedure below. The calibration consisted of four steps 
performed for each channel: (1) in-place calibration of the A/D converter and DISA 

signal conditioner; (2) static calibration for determining the King's law constants; 

(3) yaw calibration with wires in the u-w plane; and (4) recalibration to determine 
the effective angle in the measurement plane (if differenct from the u-w plane). 

Data acquisition with the calibrated system consisted of acquiring the data from the 
two channels in a buffer then dumping the buffer with identifying information to a 
floppy-disk file. 

The shorted-input reading of the A/D converter can drift a few bits from the 
nominal value of 0 and was, therefore, checked for each run. Then a reference dc 
voltage was measured with the signal conditioners bypassed. An offset value was then 
fixed on the dc offset stage of the DISA 55D26 conditioner and the offset was deduced 
by measuring the A/D value for the known reference with added offset. Finally, a 
gain was applied to the offset reference voltage through the amplifier section of the 
DISA conditioner. The effective value of the gain factor could be deduced since the 
input reference voltage and offset were accurately known. From the known calibration 
constants, the bridge output voltage could be accurately computed from the measured 
A/D converted value. 

Static calibration of the wire velocity response was performed with the wire at 
fixed orientation in a steady flow of variable velocity. Note, however, that the 
actual orientation is not yet known, but is to be determined through calibration. If 
the calibration velocity is taken to be U, then, from equation (2), 

U eff = 0 cos i)i (9) 

where \p is the angle between the calibration velocity vector and the wire. We first 
aligned the probe so that the wires lay in the u-w plane; this wire angle is called 
i); uw . With equation (1), this yields 

U = K' (E 2 - E 2 ) 1 (10) 

uw v o' 

where 

K T « K/ (cos jp ) (11) 

uw r uw 

Thus, the constants and E* were determined from a straightforward, linear, 

least-square fit of the calibration data with n specified (n is dependent on the 
calibration range; we used n = 0.45 for 5 < U < 25 m/sec). Next, ^uw was deter- 
mined via direct yaw calibration (see below); K was then computed from equation (11). 

The wires are now set at various known angles to the calibration flow in order to 
determine the "effective" wire angle \p uw . The calibration velocity was held constant 
in magnitude and direction at a value of about 70% of the maximum calibration velocity. 
If the yaw angle relative to the effective angle \p uw is denoted 6, then equa- 
tion (3) can be used to derive an equation that relates the bridge output for 6=0 
to the output for a particular value of 6; rearranging yields an expression for the 
effective wire angle \p uw : 


1 



tan \jj 


( 12 ) 


, U eff< S - 
cos(5l) , 6) 
uw sin(6 x ) 

In practice, we computed for four different values of delta of -10°, -5°, 5°, 

and 10°. These results were averaged to get the value of ^ uw * 

Now the system is calibrated for measurement in the u-w plane. When measure- 
ments in the u-v plane were desired, one further step was required. The probe stem 
was first rotated 90° to position the wires in the u-v plane. At this point, the 
effective wire angle \p uv is unknown; although \p uv would be nominally the same as 
ip uw , it can be slightly different because of slight pitching of the probe stem. 

Another static velocity calibration was performed as described above. However, this 
time n and were fixed when the King’s law was fitted to the data; K^ v was found 
by linear least-square fit, then \p uv was computed as before: 

cos ( ip ) = K/K f (13) 

v uv uv v J 

Data acquisition now took place. Identifying information regarding, for example, run 
number and probe position, was entered from the keyboard. The current flow tempera- 
ture (measured with a thermocouple and digital readout) was also entered from the 
keyboard. Then the data buffer would be filled. A few samples were used to compute 
the average bridge voltage for use in computing the temperature shift. The tempera- 
ture correction shift was written to the floppy-disk data file along with the identi- 
fying information, calibration data, and the raw data buffer; 25 kbytes were provided 
for each data file. Probe calibrations were fairly stable and repeatable for several 
hours of running, so that 60-100 data points could be reduced using the same calibra- 
tion constants with the ambient drift correction. 


Computation of Signal Statistics 

Figure 5 is a flowchart of the data-reduction algorithm for computing signal 
statistics. Starting with raw data written into floppy-disk files as described above, 
the data reduction began with construction of a look-up table from which a velocity 
could be assigned to any raw A/D voltage. This simply required that for every possi- 
ble A/D reading, the King’s law calibration be used to compute a corresponding effec- 
tive flow velocity, U e ££. Then, after adding the temperature correction shift to each 
reading, the table was entered for every raw data sample. The result for a single 
reading would be two values of U e £f — one from each channel of the crossed wires. 
Equation (6) was then used to compute the values of the instantaneous velocity compo- 
nents u and v. Once the calibration was implemented and the instantaneous-velocity 
vector components computed for each raw data point pair, the various signal statistics 
were computed. It is noteworthy that wherever reference is made to an average value, 
the average is computed using the ’’running average” formula: 

S <a(i)> = A(i - Ns) (14) 

where 

A(i) = A(i - 1) + ta(i) - A(i - l)]/i (15) 
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Average values of the various moments were computed as defined below: 

Su = S <u(i) } /Ns ^ 

Sv = S(v(i))/Ns 

Suu = S<u(i)u(i) )/Ns 

Suv = S(u(i)v(i) )/Ns V 

Svv = S<v(i)v(i) )/Ns 

Suuv = S(u(i)u(i)v(i) )/Ns 

Suw = S<u(i)v(i)v(i) )/Ns J 

Using these definitions, the signal statistics were then computed assuming nearly 
infinite sample size: 


( 16 ) 


u = Su 
v = Sv 

u' 2 = Suu - SuSu 
v' 2 = Svv - SvSv 


"\ 


y 


u'v T = Suv - SuSv 


u ,2 v* = Suuv - 2SuSuv - SvSuu + 2SvSuSu 


u ? v t2 - Suw - 2SvSuv - SuSvv + 2SuSvSv 


(17) 


J 


SAMPLE RESULTS FOR A PLANE MIXING LAYER 


Selected results of measurements made in the near-field of a plane mixing layer 
are presented in figure 7. Figure 6 depicts the situation and needed reference quan- 
tities. The mixing-layer velocity ratio was about 2:1, with a maximum velocity of 
21 m/sec. Results for both tripped and untripped initial boundary layers are pre- 
sented in figure 7. 

The profiles of mean velocity shown in figure 7(a) were fitted to the similarity 
coordinates for the developed mixing layer as recommended by, for example, Townsend 
(ref. 2). The fit results in a collapse of the mean profiles to the error function 
shape at successive streamwise locations, and the growth rate inferred fro m the 
resulting thickness parameter can be used to check the measured values of u T v f 
shown in figure 7(d). The actual values and trends of the turbulence quantities, such 
as those shown in figure 7 (b-d) , measured using the present system, compare extremely 
well with theory and data from other experiments. Full details of the measurements in 
plane mixing layers are given in reference 3. 
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5 . CONCLUDING REMARKS 


A system for rapid computerized acquisition and processing of analog signals from 
a hot-wire anemometer has been developed. Probe calibration is also implemented in 
the system. Correction for ambient temperature drift is implemented in the software. 
Complete calibration and acquisition for a typical 20-point profile requires less than 
60 min. 

Data acquired with this system in a plane mixing layer, including turbulence 
measurements up to third-order correlations, agree well with theory and existing data. 
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APPENDIX 


SOFTWARE FOR THE HP9845B DESK-TOP COMPUTER 


The HP9845B desk-top computer used included an I/O ROM and ran programs written 
in BASIC. Three programs are included: "UWIRE," a program for single hot-wire data 
acquisition and real-time data reduction; "XWIRE," for calibration, data acquisition, 
and storage of data using the crossed hot-wire CTA system; and "UVBAR," used to 
reduce the data from files written by the data acquisition program "XWIRE." 
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10 REM PROGRAM UN I RE 

20 ! PROGRAM TO ACQUIRE SINGLE-WIRE HATH USING THE LDV-R/D Cl 

SO OPTION BASE 1 

4 8 DIM T i t 1 $ [ S U ] < Y s t r f " ! i n f o r m at i o n - + r i n q s f o r d at a f i 1 e 

50 I NTEGER B1 <3000, 4 > ; dat a buffer 

6 0 I N T E G E F: S < 3 0 0 0 > , C < 2 , 3 0 @ 0 > , T < 3 000 ! * w o r d s f r o rn C I 

70 INTEGER Ns 

8 0 INTEGER; Z e- r o , E r e f , 0 f f , E q a i ri , E o f f , E z e- r o , E 1 a R 1 ay 

90 INTEGER 0bts,0ff 1 

100 INTEGER El in <20* 

110 INTEGER EOdel , Edel < 10 > 

120 INTEGER J, Jnt r 

130 INTEGER leal , I p t , N c a 1 

140 REAL F: q a i n , G a i n 

1 59 REAL V 1 < 20 > , Lie a 1 <20 > 

160 REAL N , K , E s q 

1 ? 0 R E A L U 1 c a 1 , T c a 1 , T n o w , 0 h r , A 1 p h a 

130 REAL Vhwb,Uest 

190 R E A L. U e f f < 4 0 9 6 > ! 1 o o k - u p t ab 1 e f o r U e f f cal i b r at i o n 1 2 - b i t A / D 

2 0 u REAL V v a 1 < 3 0 

2 1 0 S H 0 R T U e- f f 1 2 < 3 8 0 8 > ! d a t a a r r a y w r i 1. 1 e n t o f 1 o p f :■ y 


220 

SHORT 

Ubar 

3 0 > , U p r i < 3 0 > 

2 3 U 

PRINT 




240 

PRINTER 

IS 

0 

250 

PRINT 

l! 

■+- Hr 

■ < PROGRAM UWIRE : FULL Y -DIGITAL U-WIRE DATA 

2 6 0 

PRINT 




270 

PRINT 

n 

P R 0 1 j R A M S T R U C T 1 J R E \ 1 1 

2 8 0 

PR IN T 

ii 

1 . 

C a 1 i b r at e t h e A D e e< n y e- r t e- r •:* f t h e L D V Cl . " 

2 9 0 

PRINT 

H 

2 . 

C a 1 i b r at e t h e- r e* b e- y s . y e 1 e> c i t y . " 

3 0 0 

PRINT 

11 

d • 

C o n s t r u e t 1 o n k - ij f;. table. " 

3 1 O 

PRINT 

11 

4. 

Ac qij i re dat a and wr i t. e U t, o d i sk f i 1 e . " 

3 2 O 

PRINT 

H 

Cj 

R e p e at < 4 « > f o r e ac h d at a p o i r _ i t. t, ak e- n . !l 

3 ? 8 

PRINT 

II 

6 . 

R e d L4 c e d a t a o f f - 1 i ri e w i t h an o t. h e r* p r o g r am . 11 

3 4 u 

PR IN T 





3 5 8 ! 

5 5 0 1 * * L a 1 i b r at *= t h €• A •••'' D c h an n e 1 1 

3 70 ! 

330 PRINT 

390 PRINT "** CALIBRATION OF THE A/D CONVERTER **“ 

PRINT 

N s = 1 0 ?10s am p 1 e- s a r e ay e r aq e d at e ac h p o i n t 

hio let ,n = l 

4 >; U G 1 j b U b A d c a 1 

4 4 O E z e r o = Z e r o 

4 30 £ o f f = 0 f f 


4 6 0 

G a i n - R 

g a i n 






4 7 0 

! 







4 30 

J. Q L*i 

! * * C 

i 

a 1 i b r at e i r e y s . y e 1 o c i 

t. i,i 





5 0 0 

' 1 

. Com pi 1 e ram cal i brat } c 

n dat 

n t ab 

1 e 



510 

PRINT 







5 20 

PRINT 

CALIBRATION TO DETERMINE 

Ebr i d 

ge 

i. ,i •=. B 

Ueff 

5 3 8 

PRINT 







5 4 8 

INPUT 

"En t e r cal i b r at i o n f 1 o w 

t- e m p e 

r at u r 

t' 

i n de 

g. F:%Tcai 

55 0 

Teal =. 

5556* < leal -32 > 






560 

INPUT 

"En t e r w i r e t e m p e- r a t u r e 

res i s 

t i y i t 

i..i 

ce-eff 

i c i e n t i “ , A 1 p h a 

5 70 

INPUT 

"En t e r n o m i n a 1 o y e r h e a t 

rat i o 

u s e d 

r 

ab o u t- 

1 , 3 > : “ , Ohr 

530 

! NOTE 

i W i r e fo ar am e t e r s a r e- n e 1 

e d e d 

to do 

t. 

e m p e r*- 

at u r e e o r r e e t i o n 

5 9 8 

PRINT 

" K i n q s L aw w ill b e u s e d 

t e« 'i 

n s t u 

c t 

t. he 

look-up table - " 
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606 PRINT •' IJ = K * < E 2 “ E 0 2 > * * ( 1 •••'' N ) " 

6 10 PRINT " T h e c o n s t an t s K , E 0 "• 2 , an d N m ay b e d e t e r m i n e d f r o ru d i r e c t " 

6 2 0 P R I N T "cal i h r at i o n o r i n p u t d i r e c t 1 y . E n t e r C t o cal i b r at e , o r I “ 

6 2 O P R I N T " t o i n p u t t hi e c o n s t an t s d i r e c 1 1 y . 11 

640 INPUT " Center C or I> :",Cal$ 

658 IF Cal $ = " C " THEN GOTO 720 

66U IF Cal $= " I " THEN GOTO 680 

6 7 0 G 0 TO 64 0 

680 INPUT “Enter K : ",K 

698 INPUT " E 0 2 : ",Esq 

7O0 INPUT " N : %N 

710 GOTO 960 

? 2 O I N P U T " E n t e r n o . p o i n t s t o b e t ak e n < < = 2 8 t o t a 1 > : " , N c a 1 

73ti He = 180 1 180 samples are to be taken at each point 

740 FOR I c a 1 = 1 TO Neal 

750 PR I NT " Po i nt no . ";Ical 

7 6 0 I N P U T " E n t e r cal i b r at i o n w e 1 o c i t y : " , U c a 1 C leal > 

778 GOSUB fltod 

7 0 Eli n C leal > = 8 ! c o m p u t e ao e r ag e b r i d g e o u t p u t a 1 u e 

798 FOR Icpt=l TO Ns 

8 8 8 Eli n I c a 1 > = E 1 i n C leal > + < C C 1 , I c p t ) - E 1 i n < leal > > I e p t 

810 NEXT Icpt 

8 2 8 ! 2 . C o n e r t t o *.,« o 1 1 s 

8 3 8 V 1 C leal > = F N V b r g C E 1 i n C leal ) , E z e r o , E o f f , G a i n ) 

848 PR I NT leal ; Uc a ) C I c al > ; V 1 C I e al > 

850 NEXT leal 

868 PRINT “** CRLIBRHTION DflTR ACQUISITION COMPLETE 
87U PRINT 

8 8 8 ! 3 . P e- r f o r m K i n g •" s 1 aw f l t 

8 9 8 P PINT " P e r f o r m K i n g " s 1 aw fit — " 

980 PRINT "Channel 1 : " 

910 CALL H w c a 1 C N cal , V 1 C + > , U e a 1 C * > , N , K , E s q > 

920 PRINT 
930 ! 

948 ! * ^ C o n s t r u c t 1 o o k - u p t ac» 1 e t e> i m p 1 e m e n t t h e e a 1 i b r at i o n 

9 5 0 ! 

960 PRINT 

970 PRINT ”** LOOK-UP TABLE CONSTRUCTION AND VERIFICATION 
988 PRINT 

998 ! The matrix Ueff is a look-up table- of w a 1 u e s of velocity 

1008 ! from the King's- Law fit uersus the input w a 1 u e . 

1010 FOR . T = 1 TO 4896 

1 0 2 0 V hi w b = F N V b r g C J , E z e r o , E o f f , G a i n > 

1 8 3 8 I F V h w b 2 > E s g T H E N G 0 T 0 1 0 6 8 

1 0 4 8 U e f f < J > ~ 0 

185 8 G 0 T 0 1 0 7 8 

1 3 6 8 U e f f < J > = F N K i n g C V h w b , K , E s q , N > 

1870 NEXT J 

1030 IF C a 1 $ = *' I" THEN GOTO 1248 

1090 PRINT "Re-contruet i on of cal i brat i on data l " 

1100 PRINT 

1116 ! V e r i f y 1 o o k - u p t ab 1 e b y r e - e o n s t r u c t i n g cal i b r at i o n d at a 

1128 ! 

1 1 30 PR I NT " FT . Uae t ual El U 1 CAL " 

1148 PRINT " " 

1158 FOR Icpt=l TO Neal 
1 168 Jnt r = E 1 i n C I c pt ) 

1 1 70 U 1 c al =Uef f C Jnt r > 

1 1 88 PR I NT I c pt ; Uc al C I c pt > ; E 1 i n C I e pt > ; U leal 
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1190 
1 2 0 0 
1210 

1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1 3 1 0 
1320 
1338 
1348 
1350 
1 360 
1370 
1380 
1390 
1480 
1 4 1 0 
1420 
1430 
1440 
1450 
1 460 
14 70 
1 430 
1 4 9 0 
1 500 
1510 
1520 
1530 
1540 
1550 
1 5 6 0 
1570 
1 580 
1590 
1 6 0 0 
1610 
1 620 
1 630 
1640 
1650 
166 0 
1 670 
1 680 
1 6 9 0 
1 7 0 0 
1 7 1 0 
1720 
1730 
1 740 
1 750 
1 7 6 0 
1770 


NEXT Icpt 
PRINT 

j 

! * * Hcqui r e an d $ t o r e d at a at- s u c c e s- s i y e p o i n t s : 

! 

Isftl=0 
T n o w = T c a 1 

PRINT "** DATA ACQUISITION **" 

PRINT 

P R I N T " E n t e r r u n p ar am e- 1 e r s : " 

PRINT 

I N P U T n - N o . d at a s am p 1 e s p e r p o i n t < < 3 0 0 0 > : 11 „ N s 
RE DIM Ueffl2<Ns> 

PRINT " E n t e r r e s- p o n 5- e t o d e- t e- r rri i n e 1 1 ; .- 1 p e o f data tilt t o y r i t e ; " 
P R I N T " N - n o d at a f i 1 e w r i 1 1 e n " 

P R I N T " R - r aw d at a w r i 1 1 e n at e ac h p o i n t *' 


PRINT 

H O 

- only s urn mar 

y w r i 1 1- e- r'i at t h e e n >d o f t h e 

pro 

f i 1 e " 


INPUT 

"En 

t E‘ r N , R , c* r*- 

S : Afi 1 $ 




IF Afi 

1 $ = 

"N" THEN GOTO 

1400 




PRINT 

" E n 

t e r parent f i 

1 ename ~ for ray dat a file. 

t he 

p< o i n t 

n u m b> e r " 

PR I NT 

" i s 

appended t o 

t. h i s name . Th i s u ill bet he 

n am 

e u s e- d 

for a" 

PRINT 

"summary data fi 

1 e . " 




INPUT 

"En 

t e r f i 1 e n am e 

: " , N am e $ 




INPUT 

n _ 

E n t E- r Y i f t- 

e rri p e r a t u r e c o r - - r e c t 1 o n i s d e 

s l r e 

d : " , At 

e rri $ 


I pt = 1 

! 1 . M o e t- o n E' x t 1 o c at i o n < & e* r'i t e- r f 1 o w t E' m p it' c o r r s- c t- i n g ) 

PRINT 

PRINT "POINT NUMBER : ";lpt 

I N P U T H E n t e r V 1 o c at i o n : " , V y a 1 < I p t > 

I N PUT "E n t e r f 1 o w t e m p e- r a t u r e i ri deg. F : " , T n o w 

T n o w = . 5 5 5 6 * < T n o w - 3 2 > 

! 2 . 0 b t a i n r aw d at a 

GOSUB At od 

! 3 . E s t i m at e a*.' e r ag e b r i d g e o u t p u t y o 1 t ag e ( i f t p m p c o r r e c t i n g > 

E 1 a i, -0 

FOR 1=1 TO 10 
E 1 ay = E 1 ay + < C <1,1 > - E 1 ay > x 1 
NEXT I 

l J e =- t=Ueff (Elay > 

V 1 ay = F N V b r g < E 1 ay , Eze r o , E o f f , G a i n > 

F : R I N T "A p r o x . A x D y a 1 ue = M ; E 1 ay ; 

IF At e m $ < > " Y “ THE N G 0 T 0 1670 
R 1 ay = G a i n * E o f f + ( E 1 ay -Eze r o 
! 4 . C o m p u t e s h i f t t o cal i h r at i 

I s f t 1 = F N S f*i i f t- < A 1 fj h a , 0 h r , T c a 1 , T n o w . 

tJ e s t =• U e f f < E 1 ay + I s- f' t 1 ) 

V e st =F N V b r g < E 1 ay + I s ft 1 , Ez e r o , E o f f , 

P R I N T " S h i f t e d A x D y a 1 u e = ";E1 ay + I s f t 1 ; !i = > M ; V e s t ; " V o 1 t s 
! 5.1 m p 1 e m e- n t 1 o <:* k - u p X ab 1 e- 

Ubar < I pt > =0 

l J p r i < I p t > = 0 

FOR J=1 TO Ns 
J 1 1 t=C< 1 , J ) + I s t t 1 
U e f f 1 2 < J > = IJ e f f < J 1 1 t > 

IJ b a r < I p t > = U b ar < I p t > + < U e f f 1 2 < J > - U b« ar*- < I fj t > > x J 
IJ p r i < I p t > = U p r i < I p t > + ( U e f f ' 1 2 < J > * IJ e f f 1 2 < J > - U p r i < I p t > > x J 
NEXT J 

U p r l < I p t > = U p r i < I p t > - U b ar < I p t > * U b> ar < I p t > 

U p r i ( I fj t > = S Q R < U p r i < I f :• t > ) 


" = > H ; V 1 ay : 11 V o 1 t s 8 U = " ; U e s t 

o n f o r t e mi p £■ r at u e d r i f t 
R 1 ay ) 

G a i n ) 


14 


& 



1780 
1 790 
1 8 8 0 
1 8 1 0 
1820 
1 830 
1 840 
1850 
I860 
1870 
1 880 
1890 
1 9 0 0 
1910 
1920 
1 9 3 0 
1940 
1950 
196 0 
19 70 
1 980 

1 9 9 0 

2 0 0 0 
2 0 1 0 
2 0 2 O 
2 0 3 0 
2 0 4 0 
2 0 5 0 
2 0 6 0 
2 0 7 0 
2 0 8 0 
2 0 9 0 
2 1 0 0 
2110 
2120 
2130 
2 1 40 
2150 
2 1 60 
2 1 70 
2130 
2190 
2 2 0 0 
2210 
2 2 2 0 
2 2 3 0 
2240 
2250 
226 0 
2 2 7 0 
2 2 8 0 
229 0 
2 3 0 0 
2310 
2 3 2 0 
2330 
2340 
2 3 5 0 
2360 
2370 


P R I N T " Y = " ; V v a 1 < I p t > ; " U B A R = 11 ; U h a r < I p t > ; M U P R 1 hi E = 11 ; U p r i < I p t > 
PRINT 

IF fif i 1 *< >"R" THEN GOTO 1830 
! 6. Store data on floppy disk 

GO SUB Dfile 
PRINT 

I N P U T " E n t e r Y t o t ak e an o t h e r d at a p o i n t , el s e N : " , fl n s $ 

IF flns*="N" THEN 1900 
IF fins$<>"Y" THEN 1840 
I pt = I pt + 1 
GOTO 1450 

j 

RE DIM Ubar(Ipt) 

RED I M Upri(Ipt) 

REDIM V u a 1 (Ipt > 

IF Rfil$="S ,! THEN GOSUB Dfile 

PRINT " E n t e r o n e o f t h e f o 1 1 o w i n g t o p r o c e e d l " 

P R I N T " M ~ T ak e an o t h e r p r o f i 1 e w i t h s am e cal i b r a t 1 o n s 11 
PRINT "R - Recalibrate" 

PRINT " E - Exit program" 

I N P U T " F' r e s s M , R , o r E : " , H n s $ 

IF flns$="H" THEN GOTO 1200 
IF fi n s $ = " R " THEN GOTO 380 
IF flns$= 11 E 11 THEN GOTO 2030 
G 0 T 0 1 9 8 0 
END 

! Ml M II I I ! II I II I | I I I II I II I I I I I I ! I I I I I | I II ! I I ! I! M ill!!!! 

! ************** END OF MAIN PROGRAM UN I RE *************** 

! m M ( m j ! m i m m I? m | i j j j m m i M | m j M j i m m m i ; i ; i m j M 


D file: ! w r i t e d at a f i 1 e t o f 1 o p p y 

PRINT 

PRINT "** DATA FILE WRITE TO FLOPPY DISK ** 
PRINT 


PRINT "A t t h i s p o i n t b e s u r e t h e r •= i s a f 1 o p p y i n d r i o e 0 o f " 

P PINT " t h e 9 8 9 5 A w i t h s p ac e f o r a f i 1 e o f 5 0 , 2 5 6 - b y t e r e c o r d s . *' 


P r e s s C 0 N T t o p r o c e 

w r i t e i n f:> r o g r 


DISP "** 

PAUSE 

DISP "**** Fi 
F i 1 e T = N arciel 

I F fif i 1 *■ = " R " THEN F i 1 e$ = Narne$£, 
MASS STORAGE I S ” : HS , 0 ,0" ! 

CREATE F l 1 e $ , 5 0 ! o 


HS 


SIGN F i 1 e $ TO #1 


ed when ready 


V fi L * < Ipt > 
set floppy dri 
pen f i 1 e «-j i t h 


o e '• 9 3 9 5 A d r i v e 
160 records 256 


6 > as c 
b y t e s i- 


PRINT # 1 ; T i t 1 $ 

I F fif i 1 * = " R " THEN PR I NT # 1 ; Yu al < I pt 


PRINT #1 ; K, Esq, N 

PRINT #l;Tcal ? T ri o w , A 1 p h a , 0 h r 


PRINT # 1 ; I sf t 1 
PRINT # 1 ; N s 
IF fi f i 1 $ = " S 11 THEN MAT 
IF Af i 1 $- " S " THEN MAT 
IF Af i 1 $= " S " THEN MAT 
IF fi f i 1 $ = " R '* THEN MAT 
ASSIGN * TO #1 


PRINT # 1 ; Y a 1 
PRINT # 1 ; Ubar 
PRINT # 1 ; Upr i 
PRINT # 1 ; U e f f 1 2 
! c 1 os e d at a f i 1 e 


DISP " * * * * F i 1 e w r i t e c o rn p 1 e t e d * 
RETURN 


! ! ! ! ! ! ! ! } j j ! M ! j i | m M | m m j i i i j m m j m j m j j m | i j ;i m ; i i j ; m | j M j 


A d c a 1 > ! cal i b r a t e t h e A x D c o n v e r t e r o f t h e- L Ii V C I 

! enter the routine with I chan and Ns set 
! !! ! This s e g rn e n t i s f o r f i x i n g t h e •=• hi o r t e d - i n p u t '-, i a 1 u e 


I e f au 1 t 
:• ac hi 
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2380 PRINT 

2 3 9 0 PRINT 11 * * * Cal i b r at- e- t h e fl D c o n u e r t e r c h an n e 1 "Sic h an ; 11 * * 

2 4 0 0 B I S P " 8 h o r t i n p u t o f c h an n el " ; I c h an ; " t h e n p r e s s C 0 N T " 

2410 PfiUSE 
2420 GOSUB fltod 

2 4 3 0 ! ft e r a g e t h e 1 8 r e ad i n g 5- t o g e t t h e 2 e r o a 1 u e 
244 y Z g* f' o = y 

2450 FOR 1=1 TO Ns 

2460 Z e r o = Z e r o + L- I c h a n < I ) 

2470 NEXT I 

2 4 8 0 Z e r o = Ze r o N s 

24 9 0 P R I N T " Z e r o - 1 e v e 1 o u t p u t i s : 11 ; Z e r o ; 11 o f 12 b i t s 

2580 ! M ! Now read the reference oalue 

2518 P R I N T M N o w cal i h r at e t h e G ft I N = 1 s. DIRE C T h e f o r e p r o c e e d l n g " 

2 5 2 0 D I S P " C o n n e c t REF o 1 t ag e t o c h an " ; I c h an ; " t h e n p r e s s C 0 N T " 

2530 PAUSE 

2540 GOSUB fltod 

Z 0 5 8 E r e t~ U 

2560 FOR 1=1 TO Ns 

2578 E r e f = E r e f + C c. I c h an , I ) - E r e f ) / I 

2530 NEXT I 

2 5 9 8 P PINT " R e f e r e n c e o o 1 t ag e ap p 1 i e d t o c h an n c l ; I c h an , " i s “ J E r e f 


2 6 80 ! ! ! ! 0 f f s e t cal i b r at i o n 

2610 I* I S P " fl p p 1 y an 0 F F S E T t o c h an n e 1 " l I c h an ; M t- h e n p r e s s C 0 N T " 

2620 PAUSE 

2 6 3 0 G 0 S U B fl t o d 

264 8 0 f f = 0 

2658 FOR 1=1 TO Ns 

2 6 6 0 0 f f = 0 f f + C I c h an , I - U f f I 

2678 NEXT I 

2688 Obt s = Of f 

2 6 9 0 0 f f 1 = E r e f - 0 f f 

2 7 8 0 P PINT " 0 f f s e t *.■ a 1 u e i s : " ; 0 f f 1 

2710 ! ! ! ! S e t 1 2 fl I N n o w 

2^2 8 D I S P 11 S e t. G fl I N o n c h an n e 1 !! [ I c h an ; 11 t. P e n p r e s s L- U NT 11 

2738 PAUSE 

2748 GOSUB fltod 

2 7 5 0 E g a i n = 8 

2768 FOR 1=1 TO Ns 

2 7 7 0 E g a i n = E g a i n + (. C ( I c hi ar'i , I > -Eg ai n ) I 

2730 NEXT I 

2 7 9 0 R g a i n = < E g a l n - Z e r o > ( 0 b t s - Z e r o > 

2808 PR I NT " Gai n is : " ; Rgai n 

2818 ! !!! Reset OFFSET as desired 

2 8 2 8 D I S P ,! R e s e t 0 F F S E T o n c h an n e 1” ; I c h an ; " t h e n p r e s s C 0 N T 11 

2838 PAUSE 

2648 GOSUB fltod 

2358 Of f =0 

2868 FOR 1=1 TO Ns 

2 8 7 0 0 f f = 0 f f + < C < I c h a n , I > - 0 f f ) • - I 

2880 NEXT I 

2 3 9 0 0 f f = E g a i n - 0 f f > •■■■’ R g a i ri + ill f f 1 

2908 PR I NT " F i nal OFFSET wal ue is : 11 ; Off ; " b i t s “ 

2910 RETURN 

2928 ! i j M | | m m | | m m m m j j M m | || m i m m ! i m m ; | j m m j M | s M ! M 

2 9 3 8 fl t o d l ! S u h r o u t i n e f o r i n u t- f r - o rn t- h e- L D V - fl • Ii C I 

2 9 48 ! E n t e r r o u t i n e w i t h N s = r i o . s a rri p 1 e s 

295 8 B I S P " P r e s s C 0 N T t «:• i n i t i at e id at a ac q u s i t i o n " 

2968 PAUSE 
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RESET 10 

CONTROL MASK 10; 1 
WRITE 10 10, 5; 0 

W R I j E I 0 10,5; 1 ! s t ar t h an d s h ak e b y s e 1 1 i n g C T L 0 

N t = 4 * N s 

FOR 1=1 TO 3 
B u m m y = R E A B B I N < 1 0 ) 

NEXT I 

REDIM B 1 ( N s , 4 > 

ENTER 1 0 W F H S N t N 0 F 0 R MAT; B1 < * ) ! f as t d at a ac q u i s i t i o n 
FOR 1=1 TO Hs 
C< 1 , I > = B 1 < I , 3 >+2048 
NEXT I 

B I S P " B at a ac q u s i t i o n c o rn p 1 e t e " 

RETURN 

! ! M j m m m | m j j m ! m j | m m j j m m i m M j m m M m j M i m m m j (! j j j 

3 1 30 SUB Hue al < I NTEGER Np , REAL E ( * ) , REAL IJ f * ) , REAL N , REAL K , REAL Esq ) 

314 0 ! S u b p r o g r am t o c o m p u t e h o t - w i r e- cal i b r at i o n c o n s t an t s f o r 

3 1 5 0 ! K i n g “ s 1 aw cal i b r at i o n '•,* i a 1 i n e ar 1 e as t s q u ar e f i t 

3160 OPTION BASE 1 

3 170 I N PUT " E n t e r e x p o n e n t N < ap p r o x . 0.45) : ", N 

3130 L 8 = 0 1 initialize sums for linear 1 east -squares fit 

3190 L 2 = 0 

3 2 0 0 S 2 = 0 

3210 S 4 = 0 

3220 FOR 1=1 TO Np 

32 30 L 0 = L 0 + U < I > •■■•H 

3240 L 2 = L 2 + E < I > A 2 * U < I > -• N 

3250 S 2 = S 2 + E '• I >^2 

3 2 6 0 S 4 = S 4 + E ( I > A 4 

3270 NEXT I 

3280 B = S 2 * S 2 - N p * S 4 

3 2 90 A = < L 0 * S 2 - N p *L2 ) B 

3 3 0 0 B = < S 2 * L 2 - S 4 * L 0 > B 

3 3 1 0 K = A € 1 / N > ! s c a 1 e f ac t o r 

3 3 2 0 E s q = - B A ! e f f e c t i o e n o - f 1 o ' ...« b r i d g e o u t p u t s q u a r e d 

3 3 3 0 P R I N T " S c a 1 e f ac t o r : " ; K 

3 3 4 0 P R I NT "Ez e r o s q u ar e d : " ; E s q 

3 3 50 E r r = 0 1 c o m p u t e R M S X e r r o r o f t h e f i t 

3360 FOR 1=1 TO Np 

3 3 70 E r r = E r r + ( ( IJ < I > - K * < E ( I ) * E ( I > - E s q ) ( 1 ■ • N > > ■ U ( I > ) 2 

3330 NEXT I 

3 3 9 0 E r r = 1 0 0 * S Q R < E r r N p > 

3400 PRINT "RMS percent error of the- fit I 11 ; Err ; " " 

3410 SUEENB 

3420 ! M!!!!!! | j | j M ! M !!!! I !!!!!!! I ! I ! i !! M !!!!! I !!! : !!!!!!! M !!!!!!!! ! 

3430 BEF FNVbrg< INTEGER Ein, INTEGER Zero, INTEGER Off , REAL Gain) 

344 0 I C o rn p u t e t h e b r i d q e o u t p L4 1 u o 1 t a g e b a s e d o n cal i b r a t i o n u a 1 u e s 

3 4 5 0 R E T IJ R N < ( E i n - Z e r o ) ' G a i n + 0 f f ) * 2 8 4 0 9 6 

3460 FNENB 

3478 BEF F N Y aw < REAL Bel, REAL Ebr, REAL A, REAL N , REAL EObrs) 

3400 i m i m i t j j m i m m m ii i i i i i i i t i i m i m i t i i i i j i m t | i m t i i I m it i i M i t j i I t 

3490 BEF FNK i ng (REAL Ebr, REAL K , REAL E 0 s q , REAL N) 

3 5 0 0 ! C o m p* u t e u e 1 o c i t y b as e d o n K i n q • s L aw cal i b r at i o n c o n s t an t s 
3510 RET U R N K * < E b r 2 - E 0 s q > - • < 1 / N ) 

3520 FNENB 

3530 ! !!!!!!!(!!!!!! m !( m j } M! M !!!!!!!!!!!! I !!!!!!!!! M I !!!!!!!!!!! M ! M 

3540 BEF FNSh i ft ( REAL Alpha, REAL Ohr , RE AL T ref , REAL T, INTEGER R) 

3 5 5 0 ! C o ft i p u t e t e m p e r a t. u r e d r i f t c o r r e c t. i o n t o b r i d g e- w o 1 t a g e 


297 0 
2 9 3 0 

2 9 9 0 

3 0 0 0 
38 1 0 
3020 
3 0 3 0 
3 0 4 0 
3050 
3 O 6 0 
3070 
3 0 8 0 
3 8 9 0 
3 1 0 0 
3 1 1 0 
3120 


17 



00 CO 


3560 Cpct=Al pha*<T-Tref >/<2#<0hr-l ) > ! shift in bridge output voltage 

3570 RETURN Cpct*R 
580 FNENB 

=190 ; j ; ; ; j | ; m ; ; I;; j I!!!!!!!;; j !! i Mj j MM; j M! M I!!!!!!!!!!!!!!!!!! ! 

3600 BEF FN'Taw < REAL Be 1 , REAL Ebr,REAL FI, REAL N, REAL E0brs> 

3610 FNENB 


18 



10 

REM 

PROGRAM XWIRE2 



20 

! PROGRAM TO ACQUIRE X-WIRE DATA USING THE LD 

V-H D Cl 


3 0 

OPTION EASE 1 



40 

D I M T i 1 1 * C S 0 I , V s t r $ C 8 0 1 ! i n f o r rn at i o n s t r i n g 

s for data fi 

1 e 

50 

INTEG 

E R D 1 3 0 0 O j, 4 > ! d a t a b u f f e r 



60 

INTEG 

E R C < 2 , 3 0 0 0 > ! d a t a w o r d s f r o rn C I 



70 

INTEG 

ER Ns, I chan 



y 0 

INTEG 

E R Ze r o , Ere f , 0 f f , E g a i n , E o f f < 2) , Eze r o < 2 ) 



90 

INTEG 

E R E 1 i n < 20 > , E2 i n < 20 > 



1 0 0 

INTEG 

ER E0de 1 < 2 > , Eds 1 < 2 , 1 @ > 



110 

INTEG 

E R J n t r 9 J , N cal * I p t , leal 



120 

INTEG 

ER Rlav,R2av 



130 

INTEG 

ER Flag 



140 

REAL 

0 h r , A 1 p h a , Teal , T n o w 



150 

REAL 

F: g a i n , G a i n < 2 > 



1 60 

REAL 

V 1 ( 20 > , V 2 ( 20 > , Uc al ( 2 0 > 



170 

REAL 

N < 2 , K s t r < 2 > , K < 2 > , E s q < 2 > , C a < 2 > , A n g < 2 > 



ISO 

REAL 

N 2 < 2 > , K 2 s t r < 2 > , K 2 2 > , E 2 s q ( 2 > , C 2 a < 2 > 



190 

REAL 

V 0 d e 1 2 > , V d e 1 ( 2 > , Del ( 1 0 ) , T1 a, T2a, U leal , U 2 

cal . Tlal , T 2 a 1 


2 0 0 

PRINT 




201 

PRINTER IS 16 



2 2 8 

PRINT 

<< PROGRAM XUIRE2 : FULL Y-D I G I T AL X- 

WIFE BATA ACQ 

ISITION 

2 3 0 

PRINT 




240 

PRINT 

H PROGRAM STRUCTURE : " 



250 

PR I NT 

“ 1. Cal i br at e the H-'D converter ot the 

L D V CI. M 


260 

PRINT 

" 2 • Cali b r at e t h e K P n o b e y e . y e 1 o e i t y . 

M 


2 7 0 

PRINT 

11 3. Yaw calibration of the X probe to d 

e t c r n'i i n e w i r e 

an g 1 es 

2 8 S 

PRINT 

" 4. Acquire data and write raw data fro 

rn c h £• 1 'i: 2 t 

o d i s k 

2 9 0 

PRINT 

" 5. Repeat (5.) for each data point t ak 

en, " 


3 0 O 

PRINT 

" 6 . Red u e e d at. a o f f - 1 i n e w i t h an o t h e r p 

roar am. 4 5 * * * * * 11 


3 i 0 

PRINT 

i 




■ ■■ 2 u 
3 3 0 

| e ■* 

C a 1 i b r at e t. h e- H / D - b o t h c h an r'i e 1 s 



34 0 

! 




34 1 

PRINTER IS 0 



350 

PRINT 




3 6 0 

PRINT 

»** CALIBRATION OF THE A- D CONVERTER 



378 

PRINT 




3 8 0 

! 

1 . Channel 1 i s done f i rst 



3 9 0 

Ns =10 

! 1 0 s am p 1 e s ar e a v e r ag e d at e ae h 

p o i n t 


4 01 0 

I c h an 

= 1 



418 

GOSUB 

fldc al 



4 20 

Ezero 

( I c h an ) =Ze r o 



430 

Eoff < 

I c h an > = 0 f f 



440 

G a i n ( 

I c h an > = R g a i n 



450 

i 

2. Channel 2 is cal i brat ed next 




4 60 I Chin = 2 

4 7 0 G 0 8 U E A id c a 1 

4 8 O Eze r o I c h an .•* — Z €■ r o 

4 90 E o f f < I c h an ) = 0 f f 

5 0 0 G a i n < I c h an > = R g a i n 

510 Heal $ = " N " 

520 ! 

5 8 0 ! ** C a 1 i b r at e b o t h w i r e s at f i x e d an g 1 e y s . w e 1 o c i t y 

54 0 ! 

550 ! 1. Compile raw calibration data table 

56G PRINT 

570 PRINT "** CALIBRATION TO DETERMINE Ebridge u S . Ueff **" 
580 PRINT 


19 



5 9 S 

INPUT " Ent 

er cal i 

b r a 

it i on f 

1 o«, 

1 

t 

empe 

r 

at u r 

■ e 


i n de 

g. 


F 

■ it 

Te 

al 


600 

Teal = . 5556 

★ (Teal - 

3 2 > 


















6 1 6 

INPUT "Ent. 

e r w i r e 

t e 

m p e r at 

u r e 


r 

e s i s 

t 

i v i t 

y 


e o e f t 

i e 

•i 

e 

n t 

■ it 

3 H 

1 p h a 

620 

INPUT "Ent 

er no mi 

n a 1 

ouerh 

e at 


r 

at i o 


u £ e d 

1 

< 

about 

1 

. 

8 

) 

• it 

■ 3 

Oh 

r 

630 

! NOTE: Wi 

re para 

.met 

e r s ar 

e r* 

i e 

e 

■ded 

t 

o d o 


t 

e m p e r 

at 

u 

i r 

e 

e o r 

re 

c t i o n 

631 

PRINT " ALPHA = " 

;R1 

pha; " 

OHR 


— Ha 

” 3 

Ohr; " 



TCftL 

= 

it 


Te 

al ; 

ii 

C" 

640 

PRINT 




















650 

INPUT "Ent 

e r n o . 

p o i 

n t s t o 

be 


t 

ak e n 


( — 

2 

0 

t o t a 

1 > 


: 

ti 

j 

Nea 

1 


6 6 y 

N S = 1 0 0 

! 10 

0 •=. 

amp 1 es 

ar 

p 


t o b 

e 

t ak 

e 

n 

at e 

ae 

h 


po 

i nt 



6 ? 0 

FOR leal 

= 1 TO N 

eal 


















630 

D ISP "Po 

i n t n o . 

li a 

« 

leal 

















690 

INPUT "E 

n t e r t h 

e e 

a 1 it> r a 

t i o 

n 


e 1 o 

c 

i t v* 

: 


" , Ue a 

1 ( 

T 

e 

al 

‘j 



7 0 0 

GOSUB fit 

od 



















710 

Elm ( I 

cal >=0 


! e 

o m p 

u 

t 

e ay 

0 

rage 


b 

r i d g e 

Q 

U 

t 

pu 

t. y 

al 

ue 

720 

E 2 i n < I 

cal )=0 



















730 

FOR Ic 

pt = l TO 

Ns 


















740 

E 1 i n ( I 

cal > = E 1 

i n ( 

leal > + 

( f ! '! 

1 

j 

Icpt 

[I 

-Eli 

n 


leal > 

*;i 

I 

r 

pt. 




750 

E 2 i n( I 

cal > = E 2 

i n ( 

leal > + 

( C < 

2 

3 

I c pt 

) 

-E2 i 

n 

( 

leal > 

) y 

I 

r 

pt 




7 6 0 

NEXT I 

ept 



















770 

! 2. Co 

n y e r t t 

o o 

o 1 1 s 

















7 8 9 

VI (leal > 

= FNVbrg 

(El 

i n ( I e a 

i >, 

E 

z 

e r o ( 

1 

> , Eo 

f 

f 

1 3 U 

ai 

n 

( 

i > 

> 



7 9 9 

V2( leal > 

= F N V b r q 

( E2 

i n ( I e a 

i >, 

E 

z 

e r o ( 

2 

> Eo 

f 

f 

( 2 > , G 

ai 

n 

( 

2 - 1 

> 



8 0 0 

PRINT "P 

oint 

I c a 

l ; " u = 

ii 

u 

r 

al ( I 

r 

al > ; 



V 1 = " 

• 1 ,* 

1 

( 

Ic 

al > 

■ ii 
3 

V 2 = 

8 1 8 

NEXT lea 

1 



















820 

PRINT "** 

C R L I B R R 

TIG 

N DATA 

AC 

Q 

u 

I S I T 

I 

ON C 

0 

M 

PLETE 

* 

★ 

;l 





8 3 0 

PRINT 




















840 

! 3. Pe 

r font! K 

i ng 

s 1 aw 

f i 

t 















850 

PRINT "Per 

f o r rn K i 

n q - 

s 1 aw 

f i t 


f 

o r b 

O 

t h c 

h 

a 

n n e 1 s 

- 

- 

“ 





8 6 0 

PRINT "The 

r e s u 1 1 

S W 

i 1 1 be 

a 

S 

r 

al e 

f 

ae t o 

r 


K* = 

K •••■ 

c 

0 

i. 

R > 

an 

d" 

8 7 0 

PRINT "eff 

e c t i 1 ,* e 

n o - 

f 1 OW O 1 

u t p 

u 

t 

E 0 ■ 

2 

. R 

i 

s 

the 

e f 

f 

e 

e t 

i oe 

N.i 

i re " 

y y y 

PRINT "ang 

1 e i n t 

he 

plane < 

of 

t 

h 

e w i 

r 

ii 











8 9 8 

PR I NT 




















9 0 0 

PRINT "Cha 

n n e 1 1 

: “ 


















910 

F 1 aq=0 




















9 2 0 

IF Real $ = " 

V" THEN 

FI 

aq= 1 


I 


E 0 2 


i s f 

i 


e d f o 

r 

H 

e 

al 




930 

CALL H w c a 1 

(Neal , V 

1 <* 

5 Ucal 

( ■+■ :■ 

* 

N 

2 ( 1 > 

< 

K 2 s t 

r 

( 

1 > , E2 

- Q 

t 

1 

3 

Fla 

g> 



940 PRINT 

950 PRINT "Channel 2 


9 

6 0 

CALL H 

w c a 1 

( N 

eal 

, V2 ( ★ 

> , u 

e a 

1 ( 

★ ) 

3 N 2 

( 2 

j 

,K 

2 £ t 

r 

( 2 

> 

,E£ 

so 

( 2 > , 

FI 

ag ) 

9 

78 

! K £ t 

r i s 

K 

c 0 

S ( A > , 

W h 

er 

e 

R 

i s 

th 


e 

f f e 

c 

t i 

t. 

e w 

i r 

e an 

gl 

e 

9 

8 0 

1 w h i 

: h i 

s 

det 

e r m i n 

ed 

i n 

t. 

he 

n e 

x t 


pr 

ogr 

a 

m 

s 

e- q rn 

er 

r. 



9 

90 

« Not 

e th 

at 

th 

e e f f 

e e t 

i y 

p- 

z e 

r o - 

f 1 

o 

w 

yo 1 

t 

ag 

p 

Es 

q 

i :=• n 

€' e- 

ded 

1 

0 0 0 

! for 

t he 

i..i 

aw 

eal i b 

r at 

i o 

n « 















1 

0 1 0 

IF Rea 

1 $=" 

V" 

THEN GO 

TO 

18 

90 















1 

0 2 0 

! 






















1 

0 3 0 

! * ★ V 

aw e 

al 

i b r* 

at i on 

t, o 

d 

et 

er 

m i n 


e 

f f 

e e t. 

i 

y e 


w i r 

e 

an g 1 

P S- 


1 

040 

i 






















1 

0 5 0 

PRINT- 






















1 

0 6 0 

PRINT 


YR 

U C 

RLIBR 

RTI 

ON 

- 

P 

ROE 

E 

V 

RN 

ED 

I 

N 

N 

IRE 

P 

LANE 

* 

* " 

1 

0 7 0 

PRINT 






















1 

0 8 0 

RAD 






















1 

0 9 0 

Ns= 108 

i 

18 

0 S 

am p 1 e 

s ar e 

u 

£ e 

d i 

n 

i..i 

aw 

e a 

1 

i b 

r 

at i 

o n 




1 

1 0 0 

! 1 

. Fi 


t he 

y e 1 o 

e i t 

y 

an 

d 

get 

d 

a 

t a 

at 


de 

1 

t a= 

0 




1 

1 10 

PRINT 

"The 

n 

o m i 

n a 1 y 

aw 

an 

q 1 

p 

i s 

He 

1 

t a 

- 

t. 

h e 


e f f 

e 

t i y e 

w 

i re 

1 

120 

PRINT 

“ ang 

1 e 

i n 

the 

w i r 

p 

Pi 

an 

e i 

£ 

n 

o w 

de- 

t 

er 

rn 

i ne 

d 

by t 

h e- 

it 

1 

136 

PRINT 

" yaw 

e 

al i 

brat i 

o n . 

ii 
















1 

140 

PRINT 






















1 

158 

DISP " 

S e t 

pr 

obe 

at It 

el t 

■~L — 

8 

th 

e n 

pr 

S’ 

ss 

li: 0 

N 

T" 







1 

160 

PAUSE 























20 



1 170 

GOSUE fltod 




1180 

E 0 id e 1 ( 1 > = 0 




1 190 

E 0 d e 1 <2)=0 




1200 

FOR Icpt=l TO Ns 




1210 

E 8 d e 1 ( 1 > = E 0 d e 1 ( 1 !> + ( C ( 1 , I c p t 

> - E 0 d e 1 ( 1 

> > •••■■ I *: pt 


1320 

E 0 d *=• 1 ( 2 ) = E 0 d e 1 ( 2 ) + ( C ( 2 , I c p t 

) - E 0 d el < 2 

) > I c p t 


1230 

NEXT Icpt 



1 2 40 

V0de 1 < 1 > =FNVbrg < E0de 1 < 1 > , Ezer 

o ( 1 ) , E o f f 

< 1 ■' , G a i n ( 1 ) ) 


1250 

V 8 d e 1 < 2 > = F N V b r g < E 0 d e 1 < 2 > , E z e r 

o ( 2 ) , E o f f 

( 2 > ,Gai n ( 2 > > 


1 2 6 0 

V 0 d e 1 < 1 > = V 8 d e 1 < 1 ) 2 




1270 

V 0 d e 1 ( 2 > = V 0 d e 1 ( 2 > 2 




1271 

PRINT 




1280 

PRINT "Zero values : ch.l : " 

; V 0 d e 1 ( 1 > 

- i» - l. . 

, ‘-.rf. c. * , 

V8de 1 < 2 > ; " y 

1290 

PRINT 



1380 

! 2 . V a w t h e p r o b e t h r o u g h 

a s e r i e s 

of angles 


1 3 1 0 

Nyaw = 1 




1320 

P R I N T " N o w t. h e p r o b e w ill be 

yawed thr 

o u g h a s e r i e s 

of an q 1 es " 

1 3 30 

P PINT " af t e r w h i c h t h e ay e r ag 

e e f f e c t i 

e U i re an g 1 e 

i .= c o m p u t e d . 

1340 

PRINT 




1350 

PRINT " - Point no. " ; Nyaw 




1360 

I N P U T " E n t e r t h e w i r e an g 1 e 

Delta : " 

, Del t Nyaw) 


1370 

GOSUE fit od 




1 330 

Ede 1(1, N o an > =0 




1380 

Ede 1 ( 2 , Nyaw =0 




1 400 

FOR Icpt=l TO Ns 




1 4 1 0 

Ede 1 ( 1 , Nyaw > =Ede 1(1, Nyaw > + ( 

C <1,1 c p t ') 

-Ede 1 ( 1 , Ny aw ) > 

■Icpt 

14 2 0 

E d e 1 ( 2 , N y au > - E d e l ( 2 , N y aw ) + ( 

C < 2 , I c pt > 

-*Ede 1 ( 2 , Nyaw ? •' 

I r pt 


1430 NEXT Icpt 

1440 INPUT " R e p 1 V to do a. not her point , else- N V " , AT 
1450 IF fl$="N" THEN GOTO 1490 
14 6 0 N y a w = N y a w + 1 
1470 GOTO 1350 


1480 

! 3. 

C o 

rn p u t e 

t an q e n 

t 

ij 

f 

ay 

e 

r age 


e f f 

6‘ 

l” 

t i 

i. .i 

p 

w i 

r 

e an 

g l 

p 





1499 

PRINT 


























1 5 8 0 

PRINT " 

YflW 

CflLI 

BRRTION 

D 

ft 

TH 

s 

IJ 

MMflR 

Y 

it 















1510 

PRINT " 

PT. 

YflW 

ANGLE 



T H N ( 

fl 

1 


T RN 

( 

R2 

j 

ii 










1520 

PRINT " 

— 



__ 

- 


““ 

- 

— 

- 

— 

- 

- 


- 

“ 










1530 

T 1 a=0 


























1540 

T 2 a= @ 


























1 550 

FOR I 

c pt 

= 1 TO 

N y aw 























1560 

Vdel ( 

1 > = 

FNVbr 

g (Edel ( 

i , 

I 

c P 

t I 

j 

Ezer 

o 

( 1 > 

, 

E 

o f 

f 

( 1 

•• , 

G 

a i n ( 

1 > 

‘l 





15 70 

Vdel < 

2 ~ 

FNVbr 

g < Ede 1 < 

, 

I 

cp 

t. > 

, 

Ezer 

iii 

( 2 i 


E 

0 f 

f 

( o 

• l n 

G 

a i n < 

2 •' 






1580 

I NOTE: 

Th 

e eff 

e c t i y e 

w i 

r 

e 

an 

g 

1 e i 

s 

c 0 

m 

P 

ut 

e 

d 

at 


e ac h 

P 

o i 

nt 

-- 



1 590 

T 1 a 1 = 

FNY 

aw ( De 

1 (Icpt > 

i (i i 

% 

d 

el 

< 1 

;i 

, E2s 

q 

( 1 > 

, 

N 

2 v 

1 


V0 

d 

e 1 ( 1 

) j 






1 600 

T 1 a=T 

1 a+ 

( T 1 a 1 

-T l a :*.•••• I 

cp 

t 





















1610 

T 2 a 1 = 

FNY 

aw < D e 

1 ( I c p t > 

ii 
, *> 

d 

el 

( 2 

;i 

, E2s 

q 

r 

, 

N 

2 1 

2 

1 ? 

V 0 

d 

e 1 < 2 

’) j 






1620 

T2a=T 

2a+ 

( T 2 a 1 

- T 2 a > I 

cp 

t 





















1 630 

PRINT 

Ic 

pt ; M 

" ; Del ( 

Ic 

P 

t > 

i II 

% 


" ; t 

i 

al ; 

" 


n 


T 2 

al 









1 64 0 

NEXT 

Icp 

t 
























1650 

PRINT 


























1 6 6 0 

PRINT 

11 H 

y e r aq 

ed yalu 

es 

: 


Ta 

n 

psi 

i 

= 

it 


T 1 

a 

■ ii 


T 

art p 

s i 

2 

= 

, 

T 2 

a 

1678 

PRINT 


























1 630 

INPUT 

"R 

epl y 

C t o c h 

an 

g 

p 

t h 

e 

se , 

e 

1 se 


N 

: 


, 

Rn 

s 

$ 







1 690 

IF fln 

s$< 

\ it p M 

THEN GO 

TO 


17 

20 



















1 7 0 0 

INPUT 

11 E 

n t e r 

Tan psi 

1, 


Tan 

P 

si 2 

; 

",T 

1 

a 

J 

2 

a 










1718 

PRINT 


























1728 

C 2 a ( 1 

> = 1 

S Q R < 

1 + A B S ( T 

la 


2 

*|i 




! re 

rn 

o 

y C' 


ef 

f e 

r 

t i ye 

w 

i r 

e 

an 

g l 

e 

1 730 

K 2 ( 1 > 

— K 2 

str( 1 

> * C 2 a ( 1 

> 






! 

f r o 

m 


c a 

1 

i b 

r a 

t 

i on 

C C 1 

ns- 

t ant 

s 


1 748 

C 2 a ( 2 

> = 1 

!;! y R i 

1+flESCT 

2 a 


2 





! f o 

r 


c o 

n 

St 

r u 

C 

t i n g 

1 

o o 

k- 

up 

t 

ab 


21 



1750 

K 2 < 2 > - K 2 s t r < 2 > *C2a< 2 > 



1760 

PRINT 



1 790 

I N P U T " Y t o p e r f o r m o r t h - cal. i n p 1 an e 

of 

m e a s u r •=• rn e n t < r e p* 1 y Y o r 

, fi c a 1 $ 
1SO0 

FOR 1=1 TO 2 



1 8 1 0 

K < I > =K2 < I > 



1820 

N< I >=N2< I > 



1830 

Esq < I > = E 2 s q ( I > 



co 

© 

C a < I > = C 2 a < I > 



1350 

NEXT I 



I860 

IF Real $- u Y " THEN GOTO 640 



1870 

GOTO 1991 



1880 ! 

4 . E f f e c t i y e w i r e an g 1 e i n o r t h o g o n a 1 p 1 

an 

•=• i s c o rn p u t e d 

1890 

Ca< 1 >=K< 1 > K 2 s t r < 1 > 



1900 

C a ( 2 > = K < 2 ) •••■' K 2 s t r < 2 ) 



1 9 1 0 

PRINT 



1920 

P PINT " C a 1 ib r at i o n c o n s t an t s u s e d i n 1 o o 

k- 

u p t ab 1 e c o n s t r u t i o n : " 

1930 

P R I N T " K 1 = " ; K *: 1 > ; " K 2 = " J K < 2 > 



1 940 

PRINT 11 fil = M ;Esq(l); M R2 = ,, ;Esq<2) 



1950 

PRINT n N 1 = 11 ;N(l>; u N2 = 11 ; N < 2 > 



196 0 

p PINT 11 C 0 S 1 = " ; C a 1 ) ; " C 0 S 2 = " ; C a < 

2 ) 


1970 

R n g < 1 > = 1 8 0 * R C S < C a ( 1 > > / P I 



1980 

R n q < 2 ) = 1 8 0 * R C S < C a < 2 > ) / P I 



1 9 9 0 

PRINT M WIRE ANGLES : 1 = " ; Rng ( 1 > ; " 2 = 

n 

; Rng(2> 

1991 

PRINTER IS 16 



1 992 

PRINT " * * C R T is n o w t h e d 1 1 au 1 t p r i n t •= r 


* " 


2 0 0 0 j 

2 6 10 ! ** fi c q u i r e an d s t o r e d at a at s u c c e s s i y e p o i n t s I 

2 0 2 0 ! 


2 0 3 0 

Is ft 1 = 

0 





204 0 

Isf t2 = 

0 





2 0 5 0 

T now = T 

cal 





2 0 6 0 

PRINT 

DATA ACQUISITION 





2 0 7 0 

PRINT 






2 0 8 0 

PRINT 

" E n t e r r u n p a r a rn e t «=• r £ 

: " 




2 0 9 0 

PRINT 






2 1 0 0 

PRINT 

" - E n t €• r f i 1 e n am e f o r 

output dat 

a f i 1 e s - 

n o t t o 

exceed" 

2110 

PRINT 

" 4 c h an ac t e r s - e . g . 

DATA. " 




2120 

INPUT 

" Enter file name : ", 

Nam et 




2130 

INPUT 

" - E n t e r a 1 - 1 i n e f i 1 e 

t i t 1 €• t o r 

t hi e p r o f i 1 

e : " , T i 

t 1 t 

2 1 40 

INPUT 

" - N o . d at a s am p 1 e s p e r 

fj o i n t < < 3 

000 > : " , Ns 



2150 

INPUT 

" - E n t e r Y i f t e m p e r a t u r e c o r r e c t 

i on i s des 

i red : " 

, A t e rn $ 

2160 

Ipt = l 






2170 

! 1 

. M o y e t o n e x t 1 o c at i o n 

C enter f 

1 o w t e H i i 

f c o r r e 

c t i n >3 > 

2 1 80 

PRINT 






2190 

PR I NT 

"POINT HUMBER : ";Ipt 





2 2 0 0 

I NPUT 

" E n t e r f 1 o w t e rn p e r at u r e 

i r< deg. F : 

" , T ri o w 



2210 

T n o w = . 

5556* < T n o w - 3 2 ) 






INPUT 

" E n t e r o n e - 1 i n e s t r i n g t 

o i dent i f y 

the cur re 

n t p o i n t : " Y s t r 


2 2 3 0 ! 2 . 0 b t- a i n r a w d a t a 

2240 GOSUB fltod 

225 0 ! 3 . E s t i rn a t e a y e r a g e b r i d q e o u t. p u t. * . > o ] t, a g e < i f t e rn p c o r r e c t i n q ) 

2260 IF fltemK >"Y" THEN GOTO 2390 

2270 Ela*,' = 0 

2280 E2au=0 

2290 FOR 1=1 TO 10 

2 3 0 0 E 1 ay = E 1 ay + < C < 1 , I > - E 1 ay ) I 

2310 E 2 ay = E 2 ay + < C (2, I ) - E 2 ay > I 

2320 NEXT I 
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dr i ft 


2 3 3 0 R 1 av - G a i n *'! 1 ) * E o f t c 1 ) + ( E 1 ay - Eze r o ( 1 ) ) 

2 3 4 0 R 2 a* .' = G a i n ( 2 > * E o f f C 2 > + ( E 2 ay - E z e r* o ( 2 ) ) 

2 3 5 8 ! 4 . C o m p u t e s h i f t t o cal i b r at i o n f o r t e r 

2360 I sf t 1 =FNSh i ft < R 1 pha, Ohr ,Tcal , T now , R 1 ay > 

2 3 7 0 I s f t 2 = F N S h i f t R 1 p h a , O h r, Te al , T n o w , R 2 ay ) 

2 3 8 8 ! 5 - S t o r e d a t a o n f 1 o p p y d i s k 

2390 GOSUB Utile 
2400 PRINT 


24 1 0 

PRINT " E n t e r o n e o f t h 

e f o 1 

1 o w i n g t o p r 

oceed : " 

2420 

PRINT 11 E - 

exit the pr 

o g r am 

h 


2430 

PRINT "P - 

another dat 

a p o i 

nt , same pro 

f i 1 e name " 

2440 

PRINT "N - 

n e w p r o f i 1 e 

9 s a rn 

e cal i brat l o 

n " 

2450 

PRINT "C - 

n e w p r o f i 1 e 

5 rep 

eat cal i brat 

i o n p r o c e d u r 

246 0 

INPUT "Ente 

r E, P, N, 

o r C 

: " , flns$ 


2470 

IF flns$= ,, E‘ l 

THEN GOTO 

2550 



2430 

IF A n s $ = " P 11 

THEN GOTO 

2520 



2 4 9 0 

IF A n s $ = " N " 

THEN GOTO 

2030 



2500 

IF A n s $ = " C 11 

THEN GOTO 

350 



2510 

GOTO 2460 





2520 

I pt = I pt + 1 





2530 

2540 

GOTO 2130 

i 





255 0 
2560 

END 

i 1 i i 1 1 1 M 1 

i i i ii i i i it i 

ii m i 

1 i 1 1 1 II 1 II 1 1 

I Ii 1 1 II 1 1 1 M 


25 70 ! * * * 

2580 ! ! ! ! 

2590 D f lie: 
2 6 0 0 

2610 PRINT 
2620 PRINT 
2630 PRINT 
2640 PRINT 
2650 PRINT 
2660 PRINT 
2680 File! 
2690 PRINT 
2700 MASS 
2710 CREPT 
2720 R 3 S I G 


! M ij | m !!!!!!!!!!! M !!!!! M !!!!!!!!!! M !!!!!!!!! ! !!!!!! ! 

! * ■* * * * # •* * •+■ -+- END u F M fi I N PR U L R R M X WIRE tt*##*^##*^****# 

! | i | m j M m i m j j i m m i i ; m m j j m m M j j ; j j j m j m M j m m m j 

file: ! w r i t. e r aw X - w i r e d at a t 6 a f 1 o p p y f i 1 e f o r 1 at e r 

! r e d *.4 c t l o n < e ac h r aw R x D p a i r i s s t o r e d > 

PRINT 

PRINT + WRITE RRW DRTR FILE **“ 

PRINT 

PRINT "fit this point be sure there is a floppy in drive 0 of' 1 
PRINT "the 9 y 9 5 R with space for a fi 1e of 100 records of 256" 

P R I N T " b y t e s e ac h . P r e s s C G N T i.i h e n r e ad y t o p r o c e e d J " 

Fi 1 e $ = N am e $ & V R L. $ < Ipt > 

PRINT "File " ; F i 1 e$ ; " b e i n g w r i 1. 1 e n " 

M R S S S T O RRGE I 3 " : H 8 9 0 , 0 " ! 9 8 9 5 R f 1 o p p y d r i y e s e t as d e f au 1 t 

C R E R T E F i 1 e $ , 1 0 0 ! file i s 1 G 0 r e c o r d s o f 2 56 b y t e s 

ASSIGN File*' TO #1 


2 7 3 0 

PRINT 

# 1 ; T i t 1 * 


2740 

PRINT 

#1 ; Yst r$ 


2750 

PRINT 

# 1 : T c a 1 , T n o w , 

A 1 p h a , 0 h r 

2760 

PRINT 

#1 ; Eoff < 1 > , Ez 

e r o ( 1 ■■ , G a i n 

2 77 0 

PRINT 

#1 ; Eoff <2 > , Ez 

e r o ( 2 > j, G a i n 

2 730 

PRINT 

#1 ; K ( 1 > , Esq ( 1 

) , N( 1 > , Ca( 1 

2790 

PRINT 

# 1 ; K ( 2 > j E s q < 2 

),N(2),Ca(2 

2300 

PRINT 

#1 ; Isft 1 , Isft 


23 i 0 

PRINT 

# 1 ; Ns 



c a 1 e o n s t ant s 


2 8 2 0 M AT PRINT # 1 J C ! R aw d at a f o r b o t h w i r e s 

2330 ASSIGN * TO #1 

2 84y MASS sTuRAgE 1 5 " : T 1 5 11 ! reset t ape drive as mass storage 

2850 PRINT "** FILE WRITE COMPLETE **" 

2860 RETURN 

2378 ! ! | ! | | M | M M M M | m | m | M | i | i i m | i | | i | M m i i i i i i i m i m i i i i m i 


• e t. 

i p u t v a 1 u e 


2 8 3 0 

fide al : 

! cal i brat e 

the fi 

sD cc 

n v 

2 8 9 0 


! e n t e r t h e 

rout i 

n e w i 

t h 

2 9 0 0 
2910 

! ! ! ! 
PRINT 

T h i s s e g rn t n t- i 

s for 

f i x i 

nq 

2 9 2 0 

PRINT 

" *** Cali brat 

e the 

a . n 

c 0 1 
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2939 

2940 
2950 
2960 
2970 
2980 

2 9 9 0 
3000 

3 0 1 8 
3 0 2 0 
3 0 3 0 
3 0 4 0 
3 8 5 0 
3 0 6 0 
3070 
3080 
3 0 9 0 
3 1 0 0 
3 1 1 0 
3120 
3 1 30 
3140 
3 i 50 
3160 
3170 
3180 
3190 
3200 
3210 
3 2 2 8 
3 2 3 8 
3 24 0 
3250 
3 2 6 0 
3 2 7 0 
3 2 8 0 
3 2 9 0 
3 3 0 O 
3310 
3 3 2 0 
3 3 3 0 
3 3 4 0 
3350 
3360 
3370 
3 3 8 0 
3 3 9 0 
3400 
3410 
3420 
3430 
3 4 4 0 

34 50 
3460 
3470 
3480 
3490 
3500 

35 1 0 


H I 8 P 11 8 h o r t i n p u t o f c h an n e 1 11 ; I c h an ; " t h e n p r e s s C 0 N 
PfiUSE 

GOSUB fit od 

! H e r ag e t h e 1 0 r e ad i n g s t o g e t l h e z e r o v a 1 u e 
Z e r o = 0 

FOR 1=1 TO Ns 

Z e r o = Ze r o + < C '■ I c h an , I > - Z e r o > I 
NEXT I 

PR I NT “ Zero- 1 eve 1 out put is : 11 ; Zero ; " of 12 bits M 

! ! ! ! N o w r e a d t h e r e f e r e n •: e v a 1 u e 

B I S P ' ; C o n n e c t r e f v o 1 1 a g e t o c h an " ; I c h a n ; " t h e n p r e s 
PAUSE 

G 0 S U B A t o d 
E r e f = 0 

FOR 1=1 TO Ns 

E r e f - E r e f + ( C ■*. I c h an , I - E r e f I 
NEXT I 

F‘ PINT H R e f e r e n c e v o 1 t ag e ap pi i e d t o c hi an n el " ; I c hi an ; " 

! ! ! ! Rn offset is appl i ed and cal ibrated for channel 

D ISP 11 R p p 1 y an o f f s- e t t o c h an n e 1 “ J I c h an ; " t h e n p r e s s 
PAUSE 

GOSUB At od 
Of f = 0 

FOR 1=1 TO Ns 

0 f f = 0 f f + < C ( I c hi an , I > - 0 f f > x I 
NEXT I 

Off 1=E ref -Off 

P P I N T " 0 f f s e t v a 1 u e is : " ; 0 f t 1 

! ! ! ! A gain is calibrated - nominal values are set e 

B I S P 11 S e t G A I N o n c hi an n el 11 J I c h an ; 11 t h e n p r e s s C 0 NT" 
PAUSE 

GOSUB At od 

E g a i n = 0 

FOR 1=1 TO Ns 

E g a i n = E g a l n + C I c hi an I .> - E g a i n ) I 
NEXT I 

F' g a i n = < E g a i n - Z e r o > .*•’* < E r e f - Z e r o - 0 f f 1 > 

P R I N T " G a i r, o f c h an n el ";Ic h an ; M is " ; R g a i n 
! M ! Re s e t t hi e 0 F F S E T v a 1 u e 

B I S P 11 P e s e t t h e 0 F F S E T o n c hi an n e 1 " ; I c hi an " !l t h e n p r e s 

PAUSE 

GOSUB At od 
Of f = 0 

FOR 1=1 TO Ns 

0 f f = 0 f f + < C < I c h an , I)- 0 f f ) / I 
NEXT I 

0 f f = < E q a i n - 0 f f > • ‘ R g a i n + 0 f f 1 

PR I NT " F i n a 1 OFFSET is : H ; Off ; " b i t s " 

RETURN 

! j j i |! i i m ! | ! ! |! i m i ! ! ) ! ! j ! m i | m | M ! m | ii m i i i i M m i 

A t o d : ! S u b r o u t i n e f o r i n p u t f r o m t h e L B V - A / D C I 

! E n t e r r o u t i n e w i t hi N s = n o . s am p I e s 

B I S P " P r e s s C 0 N T t. o i n i t i at e d at a ac q u s i t i o n " 

PAUSE 
RESET 10 

CONTROL MASK 10; 1 
WRITE 10 10, 5; 0 

WRITE 10 10,5; 1 1 s t a r t h a n d s h ak e b y s e 1 1 l n g C T L 0 

Nt=4*Ns 


li 


s CUNT" 


is " ; Eref 
I chan 
C 0 NT" 


x t e r n a 1 1 y 


s CO NT" 


! 1 I It I I 
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3520 FOR 1=1 TO 3 

3 5 3 0 D u m m y = R E fl D B I N ( 1 0 > 

3540 NEXT' I 

3550 REDIM Dl(Hs,4) 

3560 ! 

3 570 E N TER 10 UFH S M t N 0 F 0 R M RT;D1 < * ) ! f as t d at a ac q u i s i t i o n 

3 5 3 0 Ii I S F 1 11 Ii at a ac q u s i t i o n c o m p 1 * t e " 

3 5 9 0 FOR 1 = 1 T 0 N s ! t r a n s f e r t h e d a t a t o i n t e g e r a r r a y s 

3 6 9 0 C < 1 , I > = D 1 < I , 3 > +2048 ! t i,i o d at a w o r d s ( L D V i s s e n d i n q 4 w o r d s t o t a 1 

36 1 0 C ( 2 , I > =D 1 ( I , 4 > +2048 ! sec ond dat a word 

3620 NEXT I 

3630 RETURN 

3640 ! !!!!!!!!!!! m !| j i {!!!!!!!!!!!!!!!!!!!!!!!!!!!!! M !!!!!!!!! f !! ! 

3650 SUE Hwc al C I N TEGER Np , PERL E < * > , RERL U < * > , RERL N , PERL K , REAL Esq , I NTEGER F 1 > 
3 6 6 0 ! S u b p r o q r am t o c o m p u t e h o t - w i r e cal i b r at i o n c o n 5- 1 an t s f o r 

3670 ! King's law calibration Mia 1 i near least square fit 

3680 OPTION ERSE 1 

3 6 9 0 I H PUT " E n t e r e x p o n e n t hi < a p fj r o x . 0.4 5 > : ", N 
3700 IF F 1 = 1 THEN GOTO 3878 

3 718 L 8 = 8 ! initial i z e s u m 5. f o r 1 i n e ar 1 e as t - s q u a r e s f i t 

3720 L 2 = 0 

3 7 3 0 S 2 = 0 

3740 S 4 = 0 

3750 FOR 1=1 TO Np 

3760 L0 = L0 + lK I > ■ ' • N 

3770 L2 = L2 + E < I > ••-2 + U < I > Ml 

3 7 8 O S 2 = S 2 + E < I j ■ 2 

3 7 9 0 S 4 = S 4 + E < I > 4 

3300 NEXT I 

3 8 1 0 D = S 2 + S 2 - N f j + S 4 

3820 R= ( L0*S£-Np*L2 > * Ei 

3 8 3 0 £ = < S 2 + L 2 - S 4 + L 0 > • " D 

3 8 4 0 K = R ( 1 7 N > ! s c a 1 e f ac t o r 

3 8 5 8 E s q = - E •••■' R let t e •: t i m e n o “ f 1 o «.-.i b r i d g e <:» u t f j u t •= q u ar e d 

3 3 6 0 G 0 T 0 3 9 4 8 
3 3 7 8 L 8 = 8 

3 8 8 8 L 1 = 0 

3890 FOR 1=1 TO Np 

3 9 00 L 8 = L 8 + U < I > + < E ( I > • ' • 2 - E s q > ( 1 / N > 

3910 L 1 = L 1 + < E < I > 2-Esq > < 2 x hi > 

3920 NEXT I 

3938 K = L 0 x L 1 

3 9 4 0 P R I N T " S c a 1 e f ac t o r K* : " ; K 

3 958 PRINT 11 Ez e r o s q u ar e d E 8 • 2 : " ; E s q 

3968 Err — 0 ! compute RMS X error of the fit- 

3970 FOR 1=1 TO Np 

3 9 8 0 E r r = E r r + < < U < I > - K + < E ( I > * E ( I >-E s q > ( 1 x N > > x U I > > - • 2 

3998 NEXT I 

4 0 0 0 E r r =100 + S Q R < E r r x N p > 

4018 PRINT "RMS percent error of the fit : " J Err j " 

4028 SUE END 

4030 i ; I I ! M < ; ; I I M I ! I ! i ! I i II i i ! i ! i ! ! ! I ! i ! i i i i it i I I I j j ; ! I t ! i ; m t I 1 1 i ! i i ! 

4848 DEF FNVbrgC INTEGER Em, INTEGER Zero, INTEGER Off, RERL Gain? 

4 0 5 8 ! C o m p u t e t h e b r i d g e o u t fj u t m o 1 t aq e h as e d o n cal i b r at i o n m a 1 u e s 

4 0 6 0 R E T U R N ( E i n - Z e r o ) x G a i n + iD f f ) + 2 8 x 4 0 9 6 

4070 FNEND 

4888 ! j !!| | j || !! | !! m ! m !! m ; t !! ii !!!!!! M !!!!!!!! M ! MM M !!!!!!!!!!!!! ! 

4098 DEF F N V aw < RERL Del, RERL Ehr f RERL A, RERL N , RERL EObrs) 

4 1 0 8 ! C o rn p u t e t h e e f f e c t i m e w i r e an g 1 e q i m e n t h e n o m i n a 1 an g 1 e D e 1 > an d 


25 


4110 ! ap p r o p r i at e c o n s t an t s f r o rn cal i b r at i o n s . y e 1 o c i t y 
4 12 0 ft n g = D e 1 * P I •• •' 1 8 0 

4 1 3 0 RET l J R N < C O S < ft n g > - < < E b r 2 - ft > / ( E © b r s -ft) ) ( 1 /■ N > > / S I N < ft n q ) 

414© FNEND 

4150 | i i m i i i m i m m I j M m i i m m | | m m j m j m i i m m i m i m i M m m i i m i i t m i i 

4160 DEF FNK i ng < REftL Ehr , REftL K , REftL E0sq , REAL N> 

4 17 0 ! C o m p u t e m e* 1 o c i t y b as e d o n K i n g s L aw cal i b r at i o n c o n s t an t s 

4 I 8 0 RET I J R N K * < E b r A 2-E0 s q > ( 1 / N ) 

4190 FNEND 

4 200 ! j j m | m | || | M | | j m ! j m i m | j i j m j | j j m j i m i m M ! M m ! i j i ; m m j j m m j i 

4218 DEF FNShi ft (REftL ft 1 pha, REftL »Dhr, REftL Tref, REftL T, INTEGER R) 

4 2 2 0 ! C o ru p u t e t e rn p e r -a t u r e d r i f t e o r r e c t i o n t o b r i d g e o 1 1 a g e 

4 2 3 9 C p c t = ft 1 p h a* < T-T r e f ) ( 2 * ( 0 h r - 1 ) ) ! '■ s h i f t i n ct r i d g e o u t fj u t m o 1 1 aq e 

4248 RETURN Cpct*R 
425© FNEND 
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16 

28 
3 8 
48 
56 
68 
70 
8 0 
96 
100 
1 18 
129 
138 

1 48 
150 
160 
178 
160 
196 

2 0 8 
2 1 8 
2 2 0 


246 
258 
2 6 0 
276 
286 


3 1 6 
3 2 0 
3 3 0 
340 
358 
3 6 0 
3 7 0 
3 8 0 

3 9 6 

4 6 8 
418 
428 
436 
446 
458 
4 6 8 
478 
4 3 6 

4 96 

5 0 6 
5 1 0 
5i6 
538 
5 40 
556 
568 
576 
530 
5 9 6 


REM PROGRAM UVBAR2 

! PROGRAM TO REDUCE X-WIRE DATA ACQUIRED WITH PROGRAM XWIRE2 
! Input is 2 channels of raw data and cal i brat i on constants 
! Output is the two components of mean velocity in the wire 
! p 1 an e as w e 1 1 as. t h e i n - p 1 an e t u r b u 1 e n c e s t r e s s e s an d t h i r d 
! order products. 

OPTION EASE 1 
DIM Tit 1 $[80] , Ystr*[803 
INTEGER Ns, Nf , Fi 1 no 
I NTEGER Eof f < 2 > , Ezero < 2 > 

INTEGER I s f t 1 , I s f t 2 
INTEGER C<2,3600> 

INTEGER I,J,Jntr 
REAL Vhwb 
REAL G a i n < 2 > 

REAL T c a 1 , T n o w , A 1 p h a , U h r 
REAL K < 2 > , Esq ( 2 > , N ( 2 > , C a < 2 > 

REAL C 1 a , C 2 a 
REAL IJeff <2, 5000 > 

S H 0 R T U e f f 1 2 (. 3 0 6 6 , 2 > 

PR I NT 

PRINT " ** PROGRAM UVBAR2 - REDUCES RAW X-WIRE DATA 
PRINT 

P R I N T 11 P r o g r am o u t 1 i n e : " 

PRINT " 1 . R e ad cal i b r at i o n an d r aw d at a f r o m s p e •: i f i e d file," 

PRINT " 2 . C o n s t r u c t 1 o o k - u p t ab 1 e an d i m p 1 e m e n t cal i b r at i o n . " 

P R I N T " 3 . C o in p u t e- r u n n i n g s u rn s f o r s t at i s t i c s u p t o t h i r d o r d e r . " 

PRINT "4. Print results." 

PRINT " NOTE: Channel 1 is assumed to be U + V wire for sign convent i on, 11 

PRINTER IS 0 

MASS STORAGE IS ":H3,0,0" 

F i 1 n o = 0 

I N PUT " E n t e r p a r e n t f i 1 e n am e < o r E t o e x i t p r o g r am ) : ",N am e t 
F i 1 &l=Name$ 

IF F i 1 e t = 11 E 11 THEN GOTO 1310 

INPUT "Enter no. of data files with parent prefix (0 if parent on 1 y ) 11 ? Nf 
IF N f = 0 THEN GOTO 418 
F i 1 no = F i 1 n o + 1 

IF F i In o > N f THEN GOTO 1290 
Fi let = N am e t i V A L t < F i 1 n o > 

ASSIGN Filet- TO #1 
READ # 1 ; T i t 1 $ 

PRINT Tit It 
READ # 1 ; V s t r t 
PRINT Yst rt 

READ # 1 ; T c a 1 , T n o w , A 1 p h a , 0 hi r 
R E A D # 1 ; E o f f ( 1 > , E z e r o ( 1 ) , G a i n ( 1 > 

R E A D # 1 J E o f f (. 2 > , E z e r o (. 2 ) , G a i n ( 2 ) 

READ # 1 ; K < 1 > , Esq < 1 > , H < 1 > , Ca< 1 ) 

READ # 1 ; K C 2 > , E s q t 2 > , N < 2 > , C a < 2 > 

READ #1 ; Isft 1 , Isf t2 
READ # 1 ; Ns 
MAT READ #1;C 

! n o w r e d u c e t h e d a t a - t h r e e s t e p s : 

! 1 . C o n s t r L4 c t & i m p 1 e m e n t 1 o o k - u p t ab 1 e 

! 2 . C o n v e r t U e f f 1,2 t o U a n d V 
! 3 . U p d a t e r u n n i n q s u m s 
IF Filr.oM THEN GOTO 721 
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6 0 0 

! 1 . C o n s t r u c t & i m p 1 e m e n t 1 o o k 

-up table 

618 

1 


620 

FOR J=1 TO 5000 


630 

FOR 1=1 TO 2 


640 

V h w b = F N V b r g < J , E z e r o < I),E o f f 

< I > , G a i n < I > > 

6 5 0 

IF V h w b 2 > E s q < I > THEN GOTO 

6 8 0 

6 6 y 

Uef f < I , J > = 0 


670 

GOTO 690 


680 

IJ e f f < I , J > = F H K i n g C V h u b , K < I > , 

Esq-:! I > , N < I > ) 

690 

NEXT I 


7 0 0 

NEXT J 


710 

DISP "LOOK-UP TABLE CONSTRUCTION COMPLETED 

720 

! I m p ] e (n e n t 1 o o k - u p 


721 

Nout 1 =0 


722 

N o u t 2 = 0 


730 

FOR J=1 TO Ns 


740 

J n t r = C ( 1 , J> + Isft 1 


741 

IF C(l, JX4095 THEN GOTO 750 


743 

Nout 1 = N o u t 1 + 1 


750 

U e f f 1 2 ( J , 1 > = U e f f < 1 , J n t r ) 


760 

J n t- r = C ( 2 , J > + I s f t 2 


761 

IF C ( 2 , J > < 4 0 9 5 THEN GOTO 770 


7 6 3 

Nout 2 = N o u t 2 + 1 


7 7 0 

IJ e f f 1 2 < J , 2 ) = U e f f < 2 , J n t r > 


730 

NEXT J 


7 9 0 

D I SP " CAL I BRAT I ON I MPLENENTAT 

ION COMPLETED 


3 0 0 ! 

3 10 ! 2 . C g n e- r t U e t f 1 , 2 t o U .an d V 


830 

3 i n 1 = 3 Q R ( 1 - C a '• 1 > ' 2 > 


840 

3 i n 2 = 3 Q R < 1 - C a < 2 > ’■ 2 > 


850 

D = C a ( 1 ) + 3 i n 2 + C a 2 > + S i n 1 


3 6 O 

A 1 = 3 i n 2 / Ii ! c h an n el 1 i s 

IJ + V wire 

8 7 0 

R2 = S i n 1 D 


8 8 0 

H 3 = C a ( 2 ) X B ! c h an n el 2 i s 

IJ - V wire- 

8 9 0 

A 4 = - C a < 1 > D 


9 0 0 

Sul =0 


910 

Sol = 0 


9 2 0 

3 u 2 = 0 


9 3 0 

3 o 2 = 0 


9 4 0 

3 u = 0 


950 

3 u 2 = 0 


9 6 0 

3 y 2 u = 0 


970 

FOR Ipt=l TO Ns 


9 3 u 

U = A 1 + !J e f f 1 2 < I pt , 1 : , +A2*Ueff 

12':. Ipt, , 2) 

9 9 0 

V = A 3 + U e f f 1 2 < Ipt , 1 > + A 4 * U e f f 

1 2 < I pt , 2 > 


1 0 0 0 ! 

1 9 1 0 ! 3 . R u n n i n g s u m s f o r s t a t i s t i c s 
1020 ! 

1 0 3 0 S u 1 = S u 1 + < U - S u 1 > I p t 

1 0 4 0 S u 2 = S u 2 + < U 2 - S u 2 > / 1 p t 

1 0 5 0 3 1 = S '. ■» 1 + < V - S 1 > I p t 

1 0 6 0 S u 2 = S o 2 + < V 2 - S o 2 > / I p t 

1 0 7 0 S u o = 3 u v + < LI * V - 3 u v > / I p t 

1 0 8 0 3 u 2 = 3 u 2 * + ( U * U * V - 3 u 2 o > 

1 0 9 0 3 o 2 u = 3 o 2 u + C U * V * V - 3 o 2 u ) 

1100 NEXT Ipt 
1 1 10 U b ar = 3 u 1 
1120 Vbar=Svl 


,ipt. 

.-•ipt 
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1130 
1 140 
1 150 
1 1 60 
1 170 
1180 
1 190 
1191 
1 192 
1193 
1200 
1220 
1240 
1260 
1270 
1280 
1290 
1 300 
1 3 1 0 
1320 
1 3 30 
1 34 0 
1350 
1 360 
i 3 7 0 
1 880 
1390 
1 490 
14 10 
1 420 
1430 
1 448 
1450 


LI p r i 2 = S u 2 - S u 1 * S u 1 

V p r l 2 = S y 2 - S y 1 * S y 1 
U y b ar = S u y ~ 5 u 1 * S y 1 

U 2 y b ar = 3 u 2 y - 2 * U b ar * S u y - V b ar * S u 2 + 2 * U b a r * U b ar * V b ar 

V 2 u b ar = S y 2 u - 2 * V b ar * S u y - IJ b ar * S y 2 + 2 * V b ar * V b a r * U t* ar 


PRINT 

PRINT " REDUCED DATA : Channel 1 taken as U+V " 
PRINT 


PRINT 

"Out -of - 

range data* 

chn 1 " ; Nout 1 ; " chn 2 11 ; N 

o u t 2 

PRINT 

"Total 5 

amp 1 e size; 

" ; N s 


PR INI- 

"HEAR 

= " ; Ubar ; " 

V ERR = " ; Vbar ; " UNITS: 

L/T" 

PRINT 

" IJPRI2 

= " ; Upr i 2; " 

VPRI2 = " ; Vpr i 2; " UNITS: 

( L ■ T > 2 " 

PRINT 

"UVBAR 

= " ; ijybar; " 

UNITS: CL T)'-2" 


PRINT 

11 U2VBAR 

= " ; U 2 y b ar ; 

" V 2 U EAR = " ; V 2 u b ar ; 11 U H I T S 

: < L T > ■*’*• 3 " 

PRINT 
GOTO 3 
INPUT 

80 

"Reply Y 

t o r e d u c e 

anot her prof i 1 e , el se N • 11 , 

An st 


IF R n s $ = M Y “ THEN 310 
MASS STORAGE IS “sTiS" 
PRINTER IS 16 
END 


j | j i ; » ! { i i 1 i M i i 1 1 M l 1 1 i 1 i M M 1 ! l i i l i i li ; i l M i l i 

i ! 1 i i 1 

i 

i i 

1 ! 

i i i i ii i i i i i f n Ti F' R n G R A M 1 1 V B A R 2 i i m i i i M ; i M i i i i 

! ; i ; ! ! 

! 

; i 

j 1 

1 1 > | t II 1 II 1 M 1 1 1 1 it 1 1 11 1 It 1 1 IE M 1 II M 1 M 1 ! 1 1 1 t M 

! i 1 1 i 1 1 

! ! ! 

ii 

! ! ! 


DEF FNKingCREAL E b r , REAL K , REAL E 0 s q , REAL N > 

! c o m p u t e y e 1 o c i t y b as e d o n K i n g '* s L aw cal i b r a i t o n 
RET U R H K * < E b r A 2-E9sq ) A < 1 / N > 

FNENH 

1 | j t | ! : ji M I | M | jl ! ! I I II M I ! M I M M M I M | i M I I I i M i M ( ! ! j M ! I I 

DEF FNVbrgC INTEGER Ein, INTEGER Zero, INTEGER Off, REAL Gain) 

RET IJ R N < ( Ei n-Ze r o > / G a i n + 0 f f > *20 4 0 9 6 
FNEHD 

i i m ! ; i m |! i m j ;; i M i i i | i m i i M i | m m m m m ;i | ; m i ;; i i ; m j i ; i ; 
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X-WIRE EQUIPMENT LIST 


ITEM 

NO. 

MFG & MODEL 

DESCRIPTION 

1 

NASA DESIGN 
(DISA 55P11) 

MINIATURE X-WIRE WITH 5mhi TUNGSTEN WIRES 
SINGLE MINIATURE HOT-WIRE WITH 5pm Pt-W SENSOR 

2 

DISA 55M10 

CONSTANT-TEMPERATURE BRIDGE 

3 

DISA 55D31 

DIGITAL VOLTMETER (AVERAGING) 

4 

DISA 55D26 

SIGNAL CONDITIONER 

5 

TEKTRONIX SC503 

10 MHz STORAGE OSCILLOSCOPE (2 CHANNEL) 

6 

TEKTRONIX PG508 

50 MHz PULSE GENERATOR 

7 

NASA DESIGN 

LDV/A-D MUX AND 4-CHANNEL A-D WITH FAST 
SAMPLE-AND-HOLD 

8 

HP98032A 

HIGH-SPEED 16-BIT PARALLEL INTERFACE 

9 

HP9845B 

DESKTOP COMPUTER WITH I/O ROM INSTALLED 

10 

HP9895A 

l 

FLOPPY DISK DRIVE 

1 



Figure 1.- Crossed-wire system hardware schematic. 
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CROSS- 


PIN ARRANGEMENT 
FOR ROTATION 
ABOUT Z-AXIS 



SPRING 
LOADED VERTICAL 
STEM FOR 
YAW MOVEMENT 



♦ x 



ELEVATION 


PLAN-VIEW 


Figure 2.- NASA crossed-wire probe and holder. 



X-WIRE 

CONDITIONED 

SIGNALS 

CHANNEL 1 — 


CHANNEL 2 — 


DIGITAL MULTIPLEXER INPUTS 16-BIT WORDS 

33 J4 J5 J6 J7 J8 


WORD 3 WORD 4 WORD 5 WORD 6 WORD 7 WORD 8 



J11 J12 J 1 3 J14 

-© Q @ © 

I 

I 


OUTPUT TO 
PDP-1 1/DR-1 1 B 


HP98032A 
' INTERFACE 


OUTPUT TO 
HP9845 98032A 


o 


POWER POWER 
FUSE INPUT 


EVENT SYNC MODE 

ii ii ii "ii "Tc^l 

A-B A C B-C A BC RANDOM 

COINCIDENCE TIME 
5 50. s 


DEAD TIME 
5-50. s 


LDV EVENT SYNCHRONIZATION , NTE RFACE CONTROL 

EVENT MODE 

L DV LDV A /D 

WORDS MULTIPLEXED 


EVENT RATE A 


COUNTER CLOCK FREQUENCY 


TRIGGER 

FROM — 4 - 
TEK PG508 


II II II IC>gf 


1 0OHi 1KHZ lOOKMz lOOKHz IMHZ 10MHz 


LDV EVENT INPUTS 

ABC 
--■© © © 


COUNTER 

OVERFLOW 

O 


□ 

□ 

□ 

□ 

EVENT RATE B 

r 

: 


□ 

EVENT RATE C 

□ 


□ 

□ 


\JX1 


— WORD STATUS — 

1 2 3 4 5 6 7 8 

OOOOOOOO 

WORD DATA 


TIME 

COUNTER 

RESET 


COINCIDENCE RATE 


012345678 9 10 11 12 13 14 15 

oooooooooooooooo 

EVENT ARM 

o 

— HANDSHAKE STATUS — 

CYC 

CTLO PFLG RDV G0 R0S1 BUSV 

o o o ltxu Jl o o o o 


in 


MANUAL EVENT RATE 
EVENT WARNING 


□ 


o 


LDV INHIBIT OUTPUTS external reset 

ABC INPUT 

© © © © 

I I 


GATE TIME 

3.1 SEC (§) i>ft 


o 

POWER 


NASA-ARC 

30- LDV COMPUTER INTERFACE MOD B 


INDICATES SELECTED OPERATING MODE 

Figure 3.- NASA LDV-A/D computer interface connections. 
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Figure 4.- Calibration and data-acquisition procedure 
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Figure 5.- Data- reduct ion algorithm 
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Figure 6.“ Schematic of plane mixing layer experimental setup 
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